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NOTEWORTHY COLLECTION 


CALIFORNIA 


ISOETES MINIMA A.A.Eaton (ISOETACEAE).— 
Modoc Co., Modoc National Forest, Warner Mts, 
Dismal Swamp area, 0.5 mi S of Oregon state line; 
41.9911°, -120.1645°; 2187 m; wet soil with short 
vegetation in a seasonal seepage area; 30 Jun 2018, 
Duncan Thomas 13332 with Wendy Madar (UC). 

Soil around the plants contains numerous mega- 
spores from last year, with long blunt spines; J. 
nuttallii A.Braun ex. Engelm. megaspores lack spines. 


Previous knowledge. Isoetes minima was described 
by Eaton (1898) from Spokane County, Washington 
collections made by Wilhelm Suksdorf. Later authors 
tended to treat it as a variety of Isoetes howellii 
Engelm., a widespread species found throughout the 
West Coast states (Pfeiffer 1922; Taylor et al. 1993). 
Based on molecular analyses, Taylor et al. (2001) re- 
elevated it to species status, while recent studies 
(Hoot et al. 2006; Larsén and Rydin 2016) revealed 
that instead of being closely related to J. howellii, I. 
minima is in a deeply divergent lineage and is most 
closely related to J. nuttallii and I. orcuttii A.A.Ea- 
ton, with which it shares the feature of three-lobed 
corms (/. howellii and relatives have bilobed corms; 
Pfeiffer 1922; Taylor et al. 1993; Freund et al. 2018). 
Isoetes minima was known only from the Suksdorf 
collections until new populations were discovered by 
A. Ceska in Washington and British Columbia 
(Ceska and Ceska 2001). 

There are no herbarium records attributed to 
Isoetes minima s.s.. 1. California in either of the 
California Consortium of Herbaria (CCH) databases 
(CCH1 and CCH 2, CCH 2019) or Consortium of 
Pacific Northwest Herbaria Specimen Database 
(CPNH 2019). There is a single specimen attributed 
to I. howellii var. minima from 1896 in CCH (2019; 
Record: SEINET10754764), yet, upon examination 
of a digital image of the specimen, none of the 
annotation labels ascribe it to /. howellii var. minima, 
only to J. howellii, so this is likely an error in the 
database. From the Washington Natural Heritage 
Program (WNHP 2005), J. minima is only found in 
two counties of Washington state and one site in 
southern British Columbia. Recent collections have 
extended its known range to western Idaho (Wash- 
ington Co.), northeastern Oregon (Wallowa, Baker, 
Umatilla, and Union cos.), southeastern Oregon 
(Warner Mountains of Lake Co.), Washington east 
of the Cascade Range (Kittitas, Okanogan, and 
Columbia cos.), with additional sites in southern 
British Columbia. However, known populations are 
widely scattered across the range and few in number 
(CPNH 2019). In both Washington and Oregon, /. 
minima is considered a taxon of concern: S1 in 


Washington (Critically imperiled; WNHP 2005; 
CPNH 2019), and Category 1 in Oregon (threat- 
ened/endangered throughout range; Oregon Biodi- 
versity Information Center 2019). 

Isoetes minima, like I. nuttallii, is a terrestrial to 
partially amphibious species. According to Eaton 
(1898) and WNHP (2005), the taxon is found in 
damp sagebrush prairie areas, and it also occurs in 
seeps and seasonally mesic-to-hydric meadows 
(CPNH 2019). This habitat is similar to that of J. 
nuttallii, but the two species are geographically 
isolated: J. minima has only been found east of the 
Cascade Range along the margins of the Columbia 
Plateau and Great basin Desert, whereas J. nuttallii is 
found west of the Cascades of Oregon and Wash- 
ington, and the Sierra Nevada and Peninsular Range 
of California. 


Significance. This collection represents a range 
extension for /. minima and is the only record from 
California (CCH 2019; CPNH 2019; Jepson Flora 
Project 2019; The University and Jepson Herbaria 
2019). A second collection was made 8.4 miles north 
of the site of this on (DT 13331, 5/30/2018, OSC 
deposition pending) in the Warner Mountains on the 
Oregon side of the state line. This is a considerable 
range expansion, as the nearest known site is near 
Baker City (245 miles) and represents an extension of 
the southern range of the species by 170 miles. 


Additional comments. There are contradicting 
accounts of the coloration of the spores in this 
species. The WNHP (2005) describe the megaspores 
of I. minima as white, whereas Thomas (2015) notes 
that they are black when wet and light gray when 
dry. The megaspores from the Dismal Swamp plants 
are off-white to medium-gray in color, but dark gray, 
nearly black when wet, yet otherwise fit the 
description provided. It is possible that the colora- 
tion of the spores is variable across the species’ range. 


—FORREST FREUND, University of California, 
Berkeley, University and Jepson Herbaria, 1001 
Valley Life Sciences Building # 2465, Berkeley, CA 
94720-2465, lentic64@berkeley.edu; DUNCAN 
THOMAS, School of Biological Sciences, Washington 
State University- Vancouver, Vancouver, WA, 98686. 
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CONSERVATION GENOMICS OF WYETHIA RETICULATA (ASTERACEAE) 


DYLAN ORION BURGE 
1887 Whitney Mesa Drive, PMB 7499, Henderson, NV g0014 
dylan.orion.burge@gmail.com 


ABSTRACT 


Wyethia reticulata Greene is a rare perennial herb found only on gabbro-derived soils of the Pine Hill 
formation in the Sierra Nevada foothills of El Dorado County, CA. Wyethia reticulata is capable of both 
vegetative (clonal) and sexual reproduction, although the former is thought to be more common than the 
latter, with sexual reproduction tied to fire. The potential dominance of vegetative reproduction has 
conservation implications for W. reticulata, as populations could appear healthy in terms of stand density, 
while losing genetic diversity due to a lack of sexual recruitment. Past genetic assessments of W. reticulata 
relied on low-resolution markers from a small number of populations and did not directly address 
conservation issues. Here, double-digest random amplified DNA (ddRAD) sequencing was used to develop a 
high-resolution genomic dataset for a sample of 168 W. reticulata samples from seven locations that represent 
every known population center of the species. The data was used to 1) calculate the frequency of asexual 
reproduction, 2) determine which populations contain unique genetic variation, and 3) test whether there are 
any ecological correlates of genetic patterns. In contrast to past genetic research, the results presented here 
indicate that clonality is not common; less than 2% of plants were identified as being derived from asexual 
reproduction. Nevertheless, some locations seem to be more affected by asexuality than others; the index of 
association varied from 4.4 to 19.4 (M = 10.9). Genetic differentiation among locations was low (F,, = 0.03- 
0.04). However, clustering identified three distinct population segments that were not previously recognized. 
Elevation was the only environmental variable significantly associated with genetic differentiation. Based on 
these results, it is recommended that conservation effort be focused on genetically isolated population 


segments. 


Key Words: California, clonal, ddRAD, DNA, endemic, genet, Pine Hill, ramet. 


Wyethia reticulata Greene (El Dorado County 
mule ears), is a rare, diploid (2n=38), perennial herb 
found in the Pine Hill region of El] Dorado County, 
California. It occurs only on soil derived from 
gabbro rock of the Pine Hill Formation (Wilson 
1986; Hunter and Horenstein 1991; Wilson et al. 
2009). Wyethia reticulata is a California Native Plant 
Society (CNPS) listed plant (California Rare Plant 
Rank 1B.2; CNPS 2019) and a US Fish and Wildlife 
Service Species of Concern (USFWS 2002). It is one 
of eight rare plants associated with gabbro-derived 
soils of the Pine Hill area, a region known for 
elevated plant endemism (Wilson et al. 2009). 

Wyethia reticulata is known for its ability to 
reproduce asexually via rhizomes (Ayres 1997; Ayres 
and Ryan 1997; USFWS 2002). In this mode of 
reproduction, each seedling can form a group of 
genetically identical individuals (hereafter referred to 
as a genet) in which all of the stems (hereafter 
referred to as ramets) are genetically identical. It is 
not known how frequently sexual versus asexual 
reproduction occurs in populations of W. reticulata, 
but mass-flowering of colonies and recruitment of 
seedlings seems to be associated with fire or 
mechanical disturbance (Ayres 1997). Nevertheless, 
recruitment is low even after disturbance (Ayres 
1997), meaning that the population biology of this 
species may be strongly affected by the long-term 
persistence of genets. This idea is supported by past 
molecular work using RAPDs combined with allo- 
zymes (Ayres and Ryan 1997, 1999). These studies 


suggest that although the analyzed populations are 
dominated by a few genets, there is high genetic 
diversity among populations (Ayres and Ryan 1997), 
even compared to a sympatric congener that 1s 
thought to have a lower rate of asexual reproduction 
(Ayres and Ryan 1999). These results are thought to 
be consistent with a species in which the persistence 
of large, asexual colonies results in the genetic 
divergence of populations. 

Although W. reticulata co-occurs with two other 
species of Wyethia (Wyethia bolanderi [A.Gray] 
W.A.Weber and W. angustifolia [DC.] Nutt.), both 
of which have the same ploidy as W. reticulata, 
neither of these species is a close relative (Moore and 
Bohs 2003). What is more, in the Pine Hill region W. 
bolanderi flowers in early spring and W. angustifolia 
in late spring, while W. reticulata does not begin to 
flower until the summer (Keil 2012), which means 
that opportunities for genetic exchange are limited. 
Although neither Wyethia species that occurs with 
W. reticulata is a close relative, three closely related 
congeners do occur in California (Moore and Bohs 
2003), including Wyethia elata H.M.Hall, the sister 
species of W. reticulata. However, these species are 
found farther south, the nearest population being at 
least 100 km south of Pine Hill (Keil 2012). Thus, it is 
unlikely that genetic exchange with other Wyethia 
species is relevant to the population genetics of the 
focal species. This idea is supported by RAPD and 
allozyme data, which failed to detect any evidence of 
genetic exchange (Ayres and Ryan 1999). 
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Throughout its limited geographic range, W. 
reticulata is under threat from development, fire 
suppression, illegal dumping, and off-highway 
vehicle use, which have led to loss of some 
populations and fragmentation of others (CNPS 
2019). As outlined above, past genetic research has 
suggested that there is a low rate of sexual 
reproduction in the species, with recruitment from 
seeds taking place mainly in the immediate after- 
math of fire or other disturbance (Ayres 1997). This 
aspect of the species’ ecology makes it especially 
vulnerable to habitat fragmentation and fire sup- 
pression; in small, isolated areas of habitat the 
populations could become dominated by just a few 
genets. In this case, a population might appear 
healthy, with high stem density in the habitat, 
despite a lack of genetic diversity. It is also known 
that W. reticulata is self-incompatible and depends 
on bees to facilitate pollen transfer (Ayers 1997). In 
isolated populations dominated by a few genets, 
sexual reproduction might fail due to a lack of 
pollen donors. Finally, variation in ecological 
conditions across the small, but topographically 
complex, Pine Hill region may have driven the 
evolution of associations between genetics and 
environment (Baack et al. 2015). Such associations 
may help populations to survive in their environ- 
ment, but could also lead to barriers to geneflow 
(Baack et al. 2015) that might curtail the spread of 
beneficial genetic variation. 

To develop and implement a conservation strategy 
for any rare species, there is a need for high- 
resolution information on genetic patterns both 
within and among populations (Allendorf et al. 
2010). Past genetic research on W. reticulata yielded 
important insights into population genetic patterns in 
this species using RAPDs and allozymes (Ayres and 
Ryan 1997, 1999). However, this work focused on 
just four populations, relied on markers with a much 
lower power of resolution than is possible with direct 
DNA sequencing, and did not directly address 
conservation concerns. In the more than 20 years 
that have elapsed since the publication of these 
studies, the advent of high-throughput DNA se- 
quencing has made it feasible to develop high- 
resolution datasets for a large number of samples 
of any species based on direct DNA sequencing 
(Mardis 2017). 

Here, restriction-associated DNA (RAD) se- 
quencing (Davey et al. 2013) is used to develop a 
large, single nucleotide polymorphism (SNP) data- 
set for all major population centers of W. reticulata 
(Fig. 1). These genetic data are used to estimate 
population genetic parameters and to determine: 1) 
the species-wide extent of clonal reproduction, 2) 
- which populations are most genetically isolated and/ 
or depauperate, and 3) whether the populations are 
adapted to local climatic conditions. Based on these 
analyses, recommendations are made for conserva- 
tion. 
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MATERIALS AND METHODS 


Sampling 


For sampling of W. reticulata, seven locations in 
the Pine Hill area were selected (Table 1, Fig. 1). 
Locations were selected to represent every major 
population center of the species and the range of 
environmental conditions under which the species is 
found. Locations were selected based on information 
obtained during preliminary visits to known popula- 
tion centers carried out between 15 and 19 June 2018. 
Preliminary visits allowed the author to select large 
populations within which to implement the sampling 
scheme outlined below. 

To prevent confusion between populations and 
subpopulations, collecting locations are referred to as 
‘locales.’ Locales are spatially discontinuous groups 
of plants with centers at least 1000 m apart and edges 
at least 100 m apart. Sampling was carried out using 
a spatially explicit, repeatable approach based on 
past research into the population biology of W. 
reticulata (Ayres and Ryan 1997) and a recent review 
of best practices for population genetic sampling of 
plants with mixed clonal and sexual reproduction 
(Arnaud-Haond et al. 2007), such as W. reticulata 

Ayres and Ryan (1997) sampled three contiguous 
patches of W. reticulata with a total area of 900 m° 
and identified a total of 9 genets in this area; this 
implies that on average, there is one genet present in 
a 100 m* patch, meaning that in an area of 
continuous cover, genets are separated, on average, 
by 5.6 m (the radius of a circle with an area of 100 
m-). Ayres and Ryan (1997) also found significant 
genetic differentiation among populations separated 
by as little as 500 m. Based on this existing 
information for the focal species, a modification of 
the stratified random sampling design recommended 
by Arnaud-Haond et al. (2007) was used, as this 
method was developed for species with patchy cover. 
A custom R script (version 3.6.2; R Development 
Core Team, R Foundation for Statistical Computing, 
Vienna, Austria) written by the author was used to 
obtain 48 random GPS points for each locale; the 
points were at least 5.6 m apart inside a 300 m 
square. During a locale visit, each of the GPS points 
was visited. At some locales, visits to the points 
required the use of protective clothing and knee pads 
in order to access plants and avoid injury from the 
dense chaparral. If a colony was present at a point, a 
sample was obtained there. This continued until all 
48 points were visited. If the 48 points yielded less 
than 24 samples, the points were re-visited in random 
order, and the nearest colony sampled, selecting a 
stem at least 5.6 m distant from the nearest sampled 
stem. This was continued until 24 samples were 
obtained. At locales other than Square Nail (Table 1, 
Fig. 1), this strategy yielded sufficient samples 
without difficulties. Due to the small number of 
plants at the Square Nail Locale, most of which were 
associated with a disused fire road (Fig. 1), an 
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Fic. 1. Sampling map for Wyethia reticulata. Grid dimensions are in degrees of latitude and longitude (WGS84 datum). 
Scale bar for locale 2344 applies to all locales other than 2342. See Table | for more information on locales. 
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TABLE 1. COLLECTING LOCALES. Code, a field collection code for each locale used as a short-hand to refer to individual 
locales; each code corresponds to an herbarium voucher (collector: D. O. Burge) deposited at the UC Davis Center for Plant 
Diversity. Latitude and longitude are given in the WGS 84 datum. Replicates, the number of biological replicates from this 
locale included in sequencing. Fire, the date of the last fire within the sampling area of the locale; ‘never’ indicates that no 
fires have been known to occur at the locale in recorded history. Soi/, indicates the soil type at the locale, from examination 
of Soil Survey Geographic (SSURGO) data for El Dorado Area, California (USDA 2019). 


Code Locale name Latitude Longitude 
2340 Pine Hill 38.7210 -120.9898 
2341 Cameron Park 38.6766 -120.9716 
2342 Square Nail 38.7726 -121.0089 
2343 Ponderosa 38.6622 -120.9524 
2344 Weber Creek 38.7617 -121.0066 
2345 Farview So./ 57 -120.9945 
2347 Tiffany Hill 38.7299 -120.9604 


additional eight stems were sampled to the south of 
the originally targeted area (Fig. 1); samples at this 
location were haphazardly selected, no less than 5.6 
m apart. 

At each locale, up to four duplicate samples were 
obtained to be used as biological replicates. Duplicate 
samples were obtained from randomly selected stems, 
and were taken from the same inflorescence on that 
stem. For all samples, ligulate petals of a single 
immature inflorescence (anthers not yet dehisced) were 
taken. Tissue was preserved in two mL of absolute 
ethanol. Flower petals were selected as the target tissue 
because preliminary research on DNA extraction 
comparing immature leaf to petal tissue (Burge 
unpublished data) suggested that only flower petals 
yielded DNA in sufficient quantities for the successful 
implementation of the ddRAD method employed here. 

Due to the highly patchy nature of W. reticulata 
populations (Ayres and Ryan 1997), in which 
demography is known to be strongly influenced by 
past fires and mechanical disturbance (Ayres and 
Ryan 1997, Ayres 1997), the resulting samples were 
not evenly distributed across the targeted areas of 
habitat (Fig. 1). What is more, the requirement for 
flowering stems occasionally meant that non-flower- 
ing stems were passed over for sampling. Both the 
location and flowering status of stems can be affected 
by both fire and mechanical disturbance history 
(Ayres and Ryan 1997; Ayres 1997). However, other 
than fire history at a locale (see below), the influence 
of non-random historical factors such as mechanical 
disturbance was not considered in the sampling due 
to the difficulty of measuring and quantifying such 
disturbances. For example, at the Square Nail locale 
few plants were found, even under chaparral that had 
re-grown since a prescribed fire that burned in 1990 
(Ayres 1997); nearly all plants (flowering or not) 
located during extensive visits to the area prior to 
sampling were located within 10 m of the fire road 
(Burge unpublished data; Fig. 1). 


Environmental Data 


Climate and soil data for each sampled plant were 
inferred for use in tests of association between 


Elevation (m) Replicates Fire Soil 
605 3 Never RgE2 
480 3 Never RgE2 
320 5 Never RgE2 
460 3 2007 RfC 
200 4 1964 RgE2 
485 3 Never RfE 
450 3 Never RgE2 


genetics and environment (see Association Analysis, 
below). The following six climatic variables were 
selected to encompass important drivers of vegeta- 
tion patterns (Flint et al. 2013): 1) climatic water 
deficit, 2) actual evapotranspiration, 3) average 
minimum winter temperature, 4) average maximum 
summer temperature, 5) annual precipitation, and 6) 
elevation. Using the R package raster (Himans 
2017), GIS data from Flint et al. (2013) were used 
to obtain the first five climatic variables for each 
sampled plant. Average minimum winter tempera- 
ture was calculated by taking the average of the 
monthly minima for December, January, and Feb- 
ruary. The average maximum summer temperature 
was calculated by taking the average of the monthly 
maxima for June, July, and August. Climatic water 
deficit and actual evapotranspiration are described 
by Hijymans (2017). Elevation was taken using a 
Garmin 12 GPS unit (Garmin Ltd., Olathe, Kansas). 

Soil data was taken from GIS data layers in the 
Soil Survey Geographic (SSURGO) database for El 
Dorado Area, California (USDA 2019). For each 
locale, the map unit name (MUNAME) for the soil 
type underlying the area was recorded. For all raw 
environmental data, see Appendix S1 (doi.org/10. 
5061 /dryad.n02v6wwt2). 

Information on past fires was obtained from the 
CAL FIRE Resource Assessment Program Fire 
Perimeters (FRAP 2019), which contains data on 
all fires of greater than ten acres in the state of 
California. The data covers all such fires since 1898. 
If a fire burned all or part of the locale, the date of 
the fire was recorded. For fires smaller than ten acres, 
records maintained by the Motherlode Field Office of 
the Bureau of Land Management were examined; no 
such fires recorded by the office had affected the 
targeted sampling areas. 


DNA Extraction and Sequencing 


Genomic DNA was extracted from the flower 
petals by Eurofins BioDiagnostics (River Falls, 
Wisconsin) using the LGC sbeadex 96 well plant 
kit. Two kits were used. In order to avoid systematic 
DNA-quality biases among locales based on extrac- 
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tion, a blocking and randomization scheme was 
employed. Samples were blocked by locale, assigning 
equal numbers of individuals from each locale to 
each of the kits, and then randomly assigning each 
individual to a well within a kit. Extraction was done 
according to the manufacturer’s instructions, with 
several proprietary modifications developed by Euro- 
fins BioDiagnostics. DNA was checked for degrada- 
tion using a 1% agarose gel, and quantified using a 
spectrophotometer. Total genomic DNA was stan- 
dardized to 30 ng/uL by Eurofins BioDiagnostics and 
sent directly to Floragenex, Inc. (Portland, Oregon) 
for further quality control, library construction, and 
DNA sequencing. 

The library preparation and sequencing methods 
employed by Floragenex followed those described by 
Lozier (2014), with the exception that the present 
work used the restriction enzyme PstI (Hipp et al. 
2014). In brief, sequence identifier barcodes and 
sequence adapters were added to the genomic DNA 
after digestion by the endonuclease PstI. Barcoded 
samples were combined and the fragments sequenced 
outwards from restriction-sites using single-end DNA 
sequencing reads 100 bp long. All DNA sequencing 
was done using a run of a HiSeq 2500 DNA 
sequencing instrument (Illumina, San Diego, CA). 
Following sequencing, DNA sequences from individ- 
ual plants were separated based on their barcodes 
using the program fastq-multx (Aronesty 2013); the 
barcodes were then removed from the sequences. 


Variant Detection 


Due to the lack of an existing reference genome for 
Wyethia or a closely related member of Asteraceae, 
detection of variation among individual plants was 
done using a ‘RAD reference’ approach (Lozier 
2014), in which a kind of reference genome was 
constructed using sequence data from the plant with 
the greatest number of unique RAD clusters. We 
employed this strategy to ensure the best reference 
genome for alignment of sequences from other 
individuals. The plant selected for use as the 
reference was 2344-22 (Table 1; Appendix S1, doi. 
org/10.5061/dryad.n02v6wwt2). In developing the 
RAD reference, custom methods developed by 
Floragenex were used to cluster identical sequences 
that had 5-500X sequencing coverage (Lozier 2014), 
which produced a preliminary assembly. The assem- 
bly was then collapsed back to separate sequences 
and these sequences were realigned against the 
preliminary RAD reference genome using the pro- 
gram BWA (Li and Durbin 2009), allowing at most 
four mismatches among reads within a cluster. The 
purpose of self-alignment was to identify and remove 
repetitive DNA regions. 

Reads for each individual plant were aligned to the 
RAD reference using BOWTIE (Langmead et al. 
2009), relying on sequence quality information to aid 
in the process of match-making, allowing a maxi- 
mum of three mis-matched bases per read, and 
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permitting alignment of each read to no more than 
one region of the reference. SAMtools (Li et al. 2009) 
was used to detect SNPs and call genotypes. Criteria 
for SNP calling required a minimum phred score of 
10, a minimum of 6X sequence coverage and a 
maximum of 25% missing data across samples. 

Following SNP calling, the data were filtered 
based on the repeatability of genotype calls, as 
revealed by the 22 biological replicates included in 
the sequencing and SNP calling. Each of the loci was 
checked for genotype concordance between repli- 
cates. A locus was removed when the biological 
replicates had different genotype calls. Such discor- 
dant loci were removed in order to ensure that only 
repeatable loci were included in the final dataset. The 
process of identifying and removing discordant loci 
was carried out using a proprietary bash script 
(Floragenex, Inc.) that employed VCFtools (Dane- 
cek et al. 2011). 


Quantification of Clonality 


Previous research indicated that W. reticulata is 
strongly clonal (Ayres and Ryan 1999), reproducing 
asexually via underground rhizomes. As is the case 
with all organisms capable of clonal reproduction, this 
can lead to one genetic individual (genet) being 
represented by multiple independent vegetative parts 
(ramets) derived from that genet (Harper 1977). Thus, 
an effort was made to quantify clonality in the dataset. 

Clonality was quantified by identifying the num- 
ber of unique multilocus genotypes (MLGs) and 
assessing whether these are the product of asexual 
reproduction. When an MLG is shared among 
multiple individuals, this does not mean that such 
individuals belong to one genet and are the product 
of asexual reproduction; the dataset may simply lack 
the power to detect the genetic differences between 
closely related plants derived from sexual reproduc- 
tion. Conversely, it is possible that multiple MLGs 
belong to one genet, with subsequent somatic 
mutation. In the latter case, the group is better 
referred to as a multilocus lineage, or MLL. Overall, 
the definition of MLLs relies on selection of a genetic 
distance threshold to group MLGs into MLLs 
(reviewed by Bailleul et al. 2016). Some researchers 
have solved this problem by collecting genetic data 
on known ramets in order to directly quantify the 
amount of genetic differentiation expected among 
ramets (Douhovnikoff and Dodd 2003). However, 
such studies are very rarely done due to the logistical 
hurdles involved. More frequently, workers have 
relied on the following categories of methodology to 
group MLG: 1) clustering algorithms (Kamvar et al. 
2015), 2) modeling asexual versus sexual reproduc- 
tion (Bailleul et al. 2016), 3) visual examination of 
frequency distributions of genetic distance (Meir- 
mans and Van Tienderen 2004; Arnaud-Haond and 
Belkhir 2007; Clark and Jasieniuk 2011), and 4) 
various combinations of these methods. In all such 
methods, however, the definition of MLLs depends, 
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TABLE 2. POPULATION GENETIC SUMMARY STATISTICS. Code, a field collection code for each locale used as a short-hand 
to refer to individual locales; each code corresponds to an herbarium voucher (collector: D. O. Burge) deposited at the UC 
Davis Center for Plant Diversity. Alleles, refers to the number of alleles inferred to segregate in each locale. Heterozygosity, 
Nei’s expected heterozygosity (Nei 1978). 74, index of association (Smith et al. 1993). rg, Index of multi-locus linkage 
disequilibrium (Agapow and Burt 2001). Latitude and longitude are in the WGS84 Datum. 


Code Locale name Alleles 
2340 Pine Hill TS 
2341 Cameron Park 7394 
2342 Square Nail Jilaldi 
2343 Ponderosa 7185 
2344 Weber Creek 7100 
2345 Farview 7093 
2347 Tiffany Hill 7475 


either directly or indirectly, on assumptions made by 
the user. 

‘Methods to quantify clonality in the present study 
were selected to minimize the number of assump- 
tions, especially with respect to population genetic or 
demographic parameters that were not directly 
measured. First, the R package poppr v. 2.8.3 
(Kamvar et al. 2014) was used to group MLGs into 
MLLs based on a-priori selection of a genetic 
distance threshold. A Euclidean distance matrix 
was created using the bitwise.dist function of poppr. 
The filter_stats algorithm was then applied to this 
matrix under default parameters, identifying three 
alternative genetic distances thresholds that could be 
used to identify MLLs. Requiring a minimum genetic 
distance to distinguish one MLL from another allows 
for small differences among MLGs, such as those 
derived from DNA sequencing error or somatic 
mutations. The m/g.filter option of poppr was used, 
under default settings, to collapse MLGs into MLLs 
based on each of the three genetic distance thresh- 
olds, farthest_neighbor (merges MLGs based on 
maximum distance between individuals in either of 
the clusters; most strict), average_neighbor (merges 
MLGs based on average distance between every pair 
of individuals between clusters; intermediate strict- 
ness), and nearest_neighbor (merges MLGs based on 
minimum distance between individuals in either 
cluster; least strict). Results of the three methods 
were similar: farthest_neighbor identified one MLL 
represented by six individuals from four locales 
(2340, 2343, 2344 & 2345); average_neighbor identi- 
fied one MLL represented by five individuals from 
three locales (2343, 2344, & 2345); nearest_neighbor 
identified two MLLs represented by four individuals 
from three locales (2343, 2344, & 2345). To test 
sensitivity to the genetic distance methods used, the 
analysis was run again using dissimilarity distance 
(diss.dist in poppr) and Net’s genetic distance (nei.dist 
in poppr). The results were identical. 

The R package RClone v. 1.0.2 (Bailleul et al. 
- 2016) was used as a graphical means to identify 
MLLs, and thereby separate ramets from genets 
among the sampled stems. The pairwise allele 
distances among the sampled stems were calculated 
(genet_dist function of RClone) and plotted as 


Heterozygosity I4 ra 
0.020 | aa 0.011 
0.022 4.37 0.003 
0.019 8.96 0.009 
0.020 B77 0.008 
0.019 19.40 0.020 
0.019 16.15 0.016 
0.022 bey | 0.004 


histograms for each local. As in Meirmans and Van 
Tienderen (2004) and Clark and Jasieniuk (2011), it 
was assumed that the frequency distributions of 
pairwise distances would each contain a wide valley 
separating stems derived from asexual reproduction 
(low end of the frequency distribution) from those 
derived from sexual reproduction (high end). How- 
ever, no such valley was identified for any locale; in 
all cases, the major cluster began at an allele distance 
of more than 1000 (Appendix S82, doi.org/10.5061/ 
dryad.n02v6wwt2). As in Meirmans and Van Tien- 
deren (2004), this cluster of greater distances was 
inferred to belong to the group of samples derived 
from sexual reproduction; additional clusters were 
present in some histograms, but these were at greater 
distances. Overall, the graphical results were taken to 
indicate that no stems were derived from clonal 
reproduction; such asexual reproduction would be 
expected to result in a separate cluster of allele 
distances at the low end of the distribution (Meir- 
mans and Van Tienderen 2004; Clark and Jasieniuk 
2011), which was not observed. 

As a final alternative method of inferring the 
degree to which plants are reproducing asexually, the 
index of association (4, Smith et al. 1993) and the 
index of multi-locus linkage disequilibrium (7g, 
Agapow and Burt 2001) were inferred for each 
population. Calculations were done in poppr (Kam- 
var et al. 2014). Both these indices express linkage 
disequilibrium. As such, both tend to be higher in 
populations where there is a higher level of asexual 
reproduction. The /,4 index is similar to rz, but the 
latter accounts for the number of loci sampled, while 
the former does not. Both of the indices suggest that 
asexual reproduction is taking place mainly in locales 
2340, 2344, and 2345 (Table 2). This supports the 
results of the MLG-based analyses above, which 
identified a small number of potential clones from 
some of the same locales. 

Although the indices of association indicate that 
asexual reproduction may be taking place in some 
locales (Table 2), the MLL results suggest that nearly 
every sampled plant is a separate genetic individual. 
Thus, it appears that asexual reproduction (clonality) 
is not common in W. reticulata, at least at the scale 
(above 5.6 m horizontal distance) that sampling was 
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done. Given the low rate of clonality (only two 
MLLs identified using the least conservative filtering 
algorithm of poppr; less than 1.2% of samples), 
asexual reproduction was ignored in subsequent 
population genetic analyses; all of the work described 
below was carried out on the complete dataset (168 
samples; Table 1), excluding the biological replicates 


Genetic Differentiation Among Locales 


Population genetic analysis was used to obtain 
genetic diversity statistics for each locale and to 
estimate the amount of differentiation among locales. 
The R packages adegenet, v. 1.3-1 (Jombart 2008; 
Jombart and Ahmed 2011) and poppr, v. 2.8.0 
(Kamvar et al. 2014, 2015) were used to calculate 
population genetic statistics. The package adegenet 
was used to calculate allelic diversity and heterozy- 
gosity for each locale. Differentiation (F,,) among 
locales was calculated using the pairwise.fst algo- 
rithm of the R package hierfstat (Goudet 2005). The 
significance of differentiation was tested using 
Goudet’s G-statistic Monte Carlo test, which was 
implemented in adegenet with 500 replicates (Jom- 
bart 2008). Isolation by distance at the level of the 
locale was assessed using partial Mantel tests carried 
out in the R package vegan v. 2.5-1 (Dixon 2003), 
controlling for the latitude and longitude of the 
locale centroid (Table 1); Pearson and Spearman 
correlation coefficients were used with 1000 replicates 
of the permutation test. As an additional means to 
quantify the relationship between genetics and spatial 
distance, the R package EcoGenetics v 1.2.1-5 was 
used to test for spatial autocorrelation; the eco.ma- 
lecot algorithm was implemented, testing for global 
autocorrelation between great circle distance and 
kinship (Nason’s Fj; Loiselle et al. 1995) 


Clustering Analysis 


The program STRUCTURE, v. 2.3.4 (Pritchard et 
al. 2000; Falush et al. 2003, 2007; Hubisz et al. 2009) 
was used to cluster the sampled individuals. 100,000 
Markov chain Monte Carlo (MCMC) replicates were 
run, discarding 20,000 replicates as burnin (Hubisz et 
al. 2009). The analysis was repeated eight times, 
varying K (the number of theoretical clusters) from 
one to eight, with 10 MCMC replicates at each K. 
The optimal K was inferred using STRUCTURE 
HARVESTER (Earl and vonHoldt 2012), according 
to the method of Evanno et al. (2005). 


Association Analysis 


Association between environment and genetics was 
tested using redundancy analysis (RDA; Forester et 
al. 2018). RDA was done using the R package vegan 
(Dixon 2003; Bourret et al. 2014). An R _ script 
developed by Forester et al. (2018) was used to impute 
missing genetic data. Unconstrained RDA was run on 
the resulting dataset to test for collinearity of the 
variables; all of the variables with a variance inflation 
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factor greater than 10 (Oksanen et al. 2018) were 
removed, which left elevation, climatic water deficit 
(hereafter ‘cwd’) and actual evapotranspiration (here- 
after ‘aet’). A constrained RDA was then carried out, 
constraining on latitude and longitude. ANOVAs 
were run to test for the significance of the full model, 
the significance of each RDA axis, and the significance 
of each of the variables. In all cases, 100 replicates of 
the permutation test were used (Oksanen et al. 2018; 
Forester et al. 2018). 


RESULTS 


Genotypes 


A total of 190 DNA samples were subjected to 
RAD sequencing, including 22 biological replicates 
(Appendix S1, doi.org/10.5061/dryad.n02v6wwt2). 
Quality was high for all samples; none were excluded 
due to low sequencing coverage or excessive missing 
data. Coverage averaged 5,536,301 reads per individ- 
ual for the 190 samples (SD = 2,030,746). Raw DNA 
sequence data is at NCBI (BioProject accession 
PRJNA549276). A total of 14,013 SNP loci were 
called (Appendix S83, doi.org/10.5061/dryad. 
n02v6wwt2). Of these, 5,930 loci passed all filters. 
No additional filters were applied; the dataset contains 
65 invariant loci (i.e., 98.9% of the loci are 
polymorphic). The filtered data is at the Dryad Digital 
Repository (Appendix S4, doi.org/10.5061/dryad. 
n02v6wwt2). 


Environmental Conditions 


Environmental conditions did not vary strongly 
across the locales (Table 1), despite the rugged nature 
of the terrain (Fig. 1). Elevation varied from 200 m at 
the Weber Creek locale, at the banks of the South 
Fork of the American River, to more than three 
times this altitude, at the summit of Pine Hill (602 m; 
Table 1). Average annual rainfall varied from 745 
mm at the Weber Creek locale, 2344, to 817 mm at 
Cameron Park, 2341 (Table 1; Appendix S1, doi.org/ 
10.5061/dryad.n02v6wwt2). The lowest climatic wa- 
ter deficit (cwd) was at Weber Creek (733 mm), and 
the highest at Square Nail (833 mm); both lowest and 
highest evapotranspiration were at Weber Creek (387 
to 476 mm), likely as a result of the varied aspect (N 
to W) at this site. 

Soils were not highly variable; all but two locales 
had the most common soil type on the Pine Hill 
formation, RgE2 (Rescue extremely stony sandy 
loam, 3 to 50 percent slopes, eroded); locale 2343 had 
RfC (Rescue very stony sandy loam, 3 to 15 percent 
slopes), and locale 2345 had RfE (Rescue very stony 
sandy loam, 30 to 50 percent slopes). 

Most of the sampled locales have never experi- 
enced wild or controlled fire, at least since records 
began to be maintained by the state of California in 
the late 19th century (Table 1). The only two 
exceptions are Ponderosa (locale 2343), which 
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experienced a fire in 2007, and Weber Creek (locale 
2344), with a fire in 1964. Although fires have burned 
in the vicinity of some of the other locales, 
potentially affecting nearby populations of the focal 
species, these fires did not burn the habitat occupied 
by the plants that were sampled in these areas. 


Population Genetics 


Allelic diversity varied from 7093 (locale 2345) to 
7475 (locale 2347; Table 2). Expected heterozygosity 
varied from 0.019 (locales 2342, 2344, and 2345) to 
0.022 (locales 2341 and 2347). Genetic differentiation 
(F,,) varied from 0.02 (locale 2340 versus 2347) to 
0.04 (locale 2343 versus 2344). Results of F,, analysis 
are at the Dryad Digital Repository (Appendix S5, 
doi.org/10.5061/dryad.n02v6wwt2). The result of 
Goudet’s G-statistic Monte Carlo test was significant 
(P = 0.001; observation: 1189.7; standard deviation 
of observation: 7.8 ; expectation: 992.7; variance: 
643.2), suggesting that there is significant hierarchical 
structuring of the locales. For tests of isolation by 
distance, the Spearman method yielded a significant 
result (P = 0.03; Mantel r statistic: 0.06), while the 
Pearson method did not (P = 0.08; Mantel r statistic: 
0.05). Spatial autocorrelation was detected by 
EcoGenetics in only three of the lowest distance 
categories: 1) 0-890 m, 2) 890-1890 m, and 3) 2770- 
3770 m (permutation test: P < 0.05; Fig. 2). 

Results of the Evanno method (Evanno et al. 2005) 
- applied to STRUCTURE output imply that three 
ancestral groups is most optimal. A summary of the 
Evanno et al. (2005) method applied to the STRUC- 
TURE output is at Dryad Digital Repository (Appen- 


dix S6, doi.org/10.5061/dryad.n0O2v6wwt2). The 
STRUCTURE plot for K = 3 (Fig. 3) shows that 
locales 2343 and 2345 are the most distinct, as they are 
dominated by two ancestral groups (Ancestral Groups 
2 and 3) that are represented in other locales, but are 
not common. The remaining locales are dominated by 
Ancestral Group 1. 


Redundancy Analysis 


RDA on environmental versus genetic data, 
controlling for latitude and longitude, yielded a 
significant overall model (F = 2.09, df = 3, P = 
0.01), but only the first axis was significant on its own 
(RDAT: Fi =i3:63,df = huPs= 00 b= RDA2: F:=51253) 
Gi fb P=OnN9 MRDAS? Fs.) 12d hoe = O25); 
Out of the three climatic variables tested, only 
elevation was significant (elev: F = 1.42, d'=1, P= 
0:04: cwd? B=19304dl = 1 P =O07aet: F = 1:302eT 
= 1, P=0.09). A plot of the first two axes (Appendix 
S7, doi.org/10.5061/dryad.n02v6wwt2) shows the as- 
sociation of climate with genetics, the primary 
differentiation occurring between the lower-elevation 
locales (2340, 2345, and 2347) and the higher elevation 
locales (2341-2344); the higher elevations locales are 
all in the vicinity of Pine Hill (Fig. 1). 


DISCUSSION 


Clonality 


The first ecological research on W. reticulata 
suggested that under natural conditions, the species 
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FIG. 3. Proportional genetic ancestry based on analysis of genotypic data in STRUCTURE. Results of analysis with the 
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reproduced sexually almost exclusively in the after- 
math of wildfire (Ayers 1997), and that in intervening 
periods, reproduction was likely asexual, via an 
underground rhizome that allows the plant to spread 
over large areas of habitat (Ayers 1997). Molecular 
genetic research (Ayers and Ryan 1997, 1999) 
supported this idea, finding that diversity was low 
within populations, but higher amongst populations. 
These results are in contrast to the results presented 
here, in which less than 2% of the 168 stems were 
identified as being derived from asexual reproduc- 
tion. These results suggest either that 1) clones are 
only present at a very fine spatial scale (less than the 
5.6 m minimum distance amongst sampled plants 
used in this study—see Materials and Methods) or 
that 2) clonal reproduction is not as common in this 
species as previously assumed. It may also be possible 
that some of the MLGs do represent ramets 
descended from one genet by asexual reproduction, 
but with sufficient somatic mutations after separation 
from the parent plant that the ramets appear to be 
separate genets. Although the methods used should 
account for such somatic mutations, these methods 
also assume a relatively low level of such change, 
while the potentially long-lived stands thought to 
occur in this species (Ayers 1997) could result in the 
accumulation of a large number of somatic muta- 
tions over the course of many years. 


Population Differentiation 


The sampled locales represent all of the major 
population centers of the species, and the complete 
range of environmental conditions under which the 
plant occurs. The results presented here suggest that 
there is population differentiation (Fig. 3) that is not 


entirely due to isolation by distance. Most signifi- 
cantly, there are two genetically distinct locales (2343 
and 2344) that are dominated by separate ancestral 
genetic groups; these locales are located at opposite 
ends of the geographic range of the species (Fig. 1), 
suggesting isolation by distance. On the other hand, 
the Square Nail locale (2342), located less than one 
km from the Weber Creek locale (2344), is not closely 
allied to plants at this locale. One potential explana- 
tion of this pattern is that the dispersal barrier posed 
by the American River Canyon has resulted in 
differentiation of these locales. A similar situation 
is found at the southern end of the species’ range, 
where the ‘Ponderosa: ‘locale"'(2343),, Pier yas 
dominated by a different ancestral group than the 
nearby Cameron Park locale (2341), despite the close 
proximity of the locations. Overall, a single ancestral 
group, Ancestral Group 1 (Fig. 3), dominates the 
core of the species’ range, centered on the higher 
elevation terrain that surrounds Pine Hill (Fig. 1). 
However, the genetic signature of the other two 
ancestral groups is significant in all of the other 
sampled locales, suggesting that genetic exchange 
among locales is more likely to take place out of, 
rather than into, the two isolated peripheral popu- 
lations (2343 and 2344) where the signature of these 
ancestral genetic groups is most common (Fig. 3). 


Ecological Correlates of Genetic Patterns 


Redundancy analysis (RDA) revealed that only 
elevation was significantly associated with genetics. 
In the RDA plot (Appendix S7, doi.org/10.5061/ 
dryad.n02v6wwt2), three of the highest elevation 
locales (2340, 2345, and 2347) are strongly differen- 
tiated from three lower elevation locales (2342-44), 
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with one intermediate elevation locale (2341) bridg- 
ing the gap. The results seem to support the findings 
of the clustering analysis (Fig. 3) in which locales 
from Pine Hill and its associated uplands (2340, 
2345, 2347; Fig. 1) are dominated by a single 
Ancestral Group, Ancestral Group 1. However, this 
same group is also at a high frequency in two other 
locales that are at lower elevations, outside the 
immediate vicinity of Pine Hill (2341 and 2342). 
These results seem to suggest that the divergent 
climate of Pine Hill and its associated uplands may 
have played a role in population differentiation in W. 
reticulata, and may continue to do so in the future. 
This may constitute a form of intraspecific local 
adaptation, of the kind that has been found to occur 
in the closely related sunflowers, Helianthus (Baack et 
al. 2015). 


Conservation Implications 


Management and conservation of rare plants 
demands knowledge of population genetic patterns 
(Allendorf et al. 2010). Data of this kind can be used 
to 1) prioritize acquisition of lands that support 
unique genetic subunits of the species (Petit et al. 
1998; Diniz-Filho et al. 2012), 2) set stringency of 
protection measures for populations depending on 
their genetic distinctiveness (Petit et al. 1998; Diniz- 
Filho et al. 2012), and 3) select germplasm for ex-situ 
conservation and re-population programs (Guerrant 
et al. 2014; Hoban and Schlarbaum 2014; Mijangos 
et al. 2015). Data presented here suggest that future 
land acquisition to support conservation of W. 
reticulata should probably prioritize populations at 
the periphery, especially at the southern and northern 
ends of the Pine Hill formation (Fig. 1), as these are 
the areas where the most genetically isolated plants 
are found (Fig. 3). These locations are also where 
populations on lands that are already under conser- 
vation should be more carefully conserved (Petit et 
al. 1998). Finally, these same populations should be a 
priority for seed-banking (Guerrant et al. 2014; 
Hoban and Schlarbaum 2014; Mijangos et al. 2015). 

As discussed above, data presented here also 
suggest that clonality may not be as common in W. 
reticulata as once thought, though some populations 
show the genetic signature of this mode of reproduc- 
tion (Table 2). This has clear conservation implica- 
tions, as the amount of clonal reproduction taking 
place determines the genetic variability of the 
populations, which in turn influences the speed at 
which populations adapt; highly clonal populations 
and species may not contain sufficient genetic 
variation to respond effectively to environmental 
changes such as habitat fragmentation, climatic 
shifts, and arrival of invasive species or pathogens 
_ (Loehle 1987; Eckert 2002). While genetic diversity is 
low in W. reticulata, the data do not support a high 
rate of asexual reproduction, at least at the scale that 
the stems were surveyed (greater than 5.6 m 
horizontal distance). However, some locales, most 
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notably Weber Creek (2344) and Farview (2345; Fig. 
1) do show a relatively strong genetic signature of 
clonal reproduction, as measured by the index of 
association (Table 2). It is not clear why these 
populations are more clonal than others; further 
research should be carried out to determine the 
conditions that predict clonality. 

Additional research is clearly needed on the 
population genetics and ecology of this rare species. 
In particular, additional data on the reproductive 
ecology of the species is needed. As past research has 
shown that sexual reproduction is strongly associated 
with fire and other kinds of major disturbance (e.g., 
bulldozing), and thus highly episodic (Ayres 1997), it 
is not clear why clonal reproduction has yet not come 
to completely dominate population genetic patterns 
in the species. Fine-scale population genetics (below 
5.6 m horizontal distance) on continuous stands of 
the species would help to unravel the scale at which 
clonality might occur. It would be especially inter- 
esting to excavate underground rhizomes to confirm 
which stems are ramets, and to compare the genetic 
profiles of these stems to the seeds that they produce, 
as this would help to place a threshold on the 
expected amount of genetic divergence among stems 
derived from asexual versus sexual reproduction 
(Arnaud-Haond et al. 2007a). More experiments on 
demographic response to fire and other disturbance 
should also be carried out, as the results of such 
research would lead to the development of manage- 
ment tools, for example disturbance to promote 
sexual reproduction, that could be used to maintain 
genetic diversity in this rare species. 
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ABSTRACT 


Oak leaf morphology is highly variable. Leaf shape in oaks has been shown to covary with climate, and leaf 
dissection is negatively correlated with hydraulic resistance. Such patterns may reflect local adaptation of 
genotypes or plastic responses to local environments. We tested two hypotheses regarding variation in leaf 
shape in Quercus lobata Née, a widely distributed oak tree endemic to California, using subsamples of trees 
from a large provenance trial involving source populations throughout the species range. First, we tested the 
hypothesis that variation in leaf dissection (perimeter - area ') and specific leaf area could be explained by 
adaptation to climate. If so, we predicted that the progeny of trees originating from contrasting climates 
would also contrast in leaf dissection or specific leaf area when grown in a common environment. Second, we 
tested the hypothesis that variation in leaf dissection is associated with photosynthetic rate. Because dissection 
affects hydraulic resistance and heat transfer, more dissected leaves are thought to experience lower water 
stress. If so, we predicted that individuals with more dissected leaves would maintain higher rates of 
photosynthesis than those with less dissected leaves, especially during hot and dry conditions. We found no 
association of leaf traits with maternal climate, but we observed a positive relationship between leaf dissection 
and photosynthetic rates. This result suggests variation in leaf shape may have functional consequences and 
influence how valley oaks cope with environmental stress in California’s Mediterranean climate. 


Key Words: adaptation, climate, common garden, dissection, leaf shape, photosynthesis, Quercus lobata, 


water stress. 


Leaves are the primary organs of carbon assimi- 
lation, water relations, and energy balance in plants, 
and thus leaf traits that affect fitness, such as 
morphology, may be under strong selection, with 
optimal leaf morphology depending on local envi- 
ronmental conditions (Givnish 1978; Wright et al. 
2005; Ferris 2019). The shape of leaves can affect 
plant thermoregulation, vulnerability to drought, 
resistance to herbivory, and light interception, 
among other processes (Nicotra et al. 2011; Ferris 
2019), any of which may have consequences for plant 
survival and reproduction. The potential for local 
adaptation in leaf traits should be especially pro- 
nounced in widely distributed species that experience 
substantial variation in local climatic conditions 
(Valladares et al. 2000; Royer 2012; Ramirez- 
Valiente et al. 2017a). 

Oaks (genus Quercus) are well known for harbor- 
ing substantial intraspecific variation in leaf mor- 
phology, and this variation may reflect different 
ecological strategies (Jensen 1990; Valladares et al. 
2000). Previous studies, for example, have shown 
that specific leaf area and leaf shape of some oak 
species covary with climatic variables, including 


temperature and precipitation (Kolb et al. 1990; 
Valladares et al. 2000; MacDonald 2017) and 
elevation (Albarran-Lara et al. 2015). Ramirez- 
Valiente et al. (2017b), for example, reported that 
genetically based variation in functional leaf traits 
among populations of Q. oleoides Schltdl. & Cham., 
including specific leaf area and lamina size, resulted 
from adaptation to precipitation regimes. In addi- 
tion, Siso et al. (2001) found that leaf dissection is 
negatively associated with hydraulic resistance in 
several North American Quercus species, indicating a 
functional consequence of leaf dissection that may 
drive adaptation. Leaf dissection appears to improve 
heat transfer (Gurevitch and Schuepp 1990), which 
likely is adaptive for plants exposed to high 
temperatures. Taken together, studies of oaks 
suggest that leaf morphological variation adaptively 
mediates responses to the environment (e.g., Abrams 
1990; Cavender-Bares et al. 2004; Riordan et al. 
2016). 

In this study, we tested two hypotheses regarding 
variation in leaf morphology of Quercus lobata Nee, 
California’s valley oak, in the context of a large 
common garden experiment (Delfino-Mix et al. 
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2015). Leaf morphology in valley oak has been found 
to vary with elevation (Albarran-Lara et al. 2015) 
and to differ genetically among populations (Mac- 
Donald 2017). In addition, the broad topographic 
and geographic distribution of valley oak, including 
the Central Valley, Sierra Nevada foothills, and 
Costal Ranges of California (Tyler et al. 2006) 
suggests that the substantial variation in climate 
across its range may have imposed selection on 
important morphological and physiological traits 
(MacDonald 2017; Browne et al. 2019). First, we 
tested the idea that variation in leaf morphology is 
explained by adaptation to the climate in which a 
population evolves. If so, we predicted that trees 
grown from acorns originating in contrasting cli- 
mates would differ in leaf morphology even when 
grown in a common garden. Second, we tested the 
hypothesis that variation in leaf morphology is 
associated with photosynthetic rate. Higher photo- 
synthetic rates are positively associated with individ- 
ual fecundity (e.g., Arntz et al. 2000) and population 
growth (e.g., Westerband and Horvitz 2017) in many 
plants (reviewed in Ackerley et al. 2000). Siso et al. 
(2001) found that more dissected leaves of several 
Quercus species experienced reduced hydraulic resis- 
tance, which may allow individuals with dissected 
leaves to have improved water balance, especially 
during periods of drought stress. If leaf dissection 
ameliorates water stress, we predict that more 
dissected leaves would have higher carbon assimila- 
tion rates (a measure of photosynthesis), particularly 
during periods of water stress, given the demands of 
photosynthesis for water use. 


METHODS 


Study Species 


Valley oak (Quercus lobata |[Fagaceae]) is an 
iconic, keystone species within California, and it 
provides a useful study system because of its broad 
geographic range, as well as the availability of 
extensive background information. Valley oak pop- 
ulations harbor fine-scale genetic structure (Dutech 
et al. 2005), leading to the potential for local 
adaptation to environmental conditions. This struc- 
ture is likely a consequence of relatively restricted 
effective distances of pollen-mediated and seed- 
mediated gene flow (Sork et al. 2002; Grivet et al. 
QOS Griver “et lal." 20098" Sork— etal 2095)? Tn 
addition, genetic structure for putatively neutral loci 
(Sork et al. 2010; Gugger et al. 2013), as well as non- 
neutral loci (Sork et al. 2016) has been shown to be 
correlated with climate gradients, which supports the 
idea that this species has adapted to climate (see also 
Browne et al. 2019). Moreover, valley oak leaves 
- appear to harbor substantial morphological varia- 
tion, from genetic as well as environmental sources 
(MacDonald 2017). A range-wide field study of 
valley oaks reported geographic variation in mor- 
phological as well as physiological traits (R. Coria, 
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FIG. 1. Map of California showing the localities from 
which Quercus lobata acorns were collected (open hexagons) 
to establish the two provenance trials (open stars) in Chico 
and Placerville. Localities of family lines subsampled from 
high (triangle), median (square), and low (circle) climate 
categories indicated with shaded symbols (see Methods). 


K. D. Gaddis, J. M. Espelta, and V. L. Sork, 
unpublished data). Valley oak is currently threatened 
throughout its range (Tyler et al. 2006; Mclaughlin 
and Zavaleta 2012) due to extensive habitat loss 
(Whipple et al. 2011), and the species is predicted to 
experience a substantial reduction in available 
habitat within the next century due to climate change 
(Kueppers et al. 2005; Sork et al. 2010). In addition, 
valley oak is sensitive to water stress, especially as 
seedlings (Meyer 2002). 


Experimental Design 


This study was conducted using a large fully 
replicated provenance test established at two loca- 
tions, in which acorns collected from maternal trees 
were grown in a common environment (for details of 
study see Delfino-Mix et al. 2015). In fall 2012, 20 
acorns were collected from each of 672 mature trees 
at 95 localities (provenances) across the native range 
of valley oak. On 14 and 15 November 2012, up to 16 
surface-sterilized acorns within a maternal family 
(total 10,326) were planted into containers and 
reared in a lath house under 50% shade cover. In 
late 2014 and early 2015, five saplings (maternal 
families) from each sampled mature tree were 
transplanted into large field plots using a randomized 
block design at each of two locations: (1) the Chico 
Seed Orchard (“Chico”) in Chico, CA and (2) the 
Institute of Forest Genetics (“Placerville”) in Placer- 
ville, CA (Fig. 1). The two sites are quite distinct 
climatically, but both fall within the range of climate 
variation across the sites from which acorns were 
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collected (Delfino-Mix et al. 2015). Chico is a 
warmer, drier climate (68 m elevation; Mean Annual 
Temperature [MAT] 16.4 °C; Mean Annual Precip- 
itation [MAP] 64.1cm), while Placerville has a cooler, 
wetter climate (838 m elevation; MAT 14.5°C; MAP 
103.8cm; Wang et al. 2016). Both tests were irrigated 
periodically during summers during our study, and 
competing vegetation was controlled with mowing 
and herbicide. See Delfino-Mix et al. (2015) for more 
information about sampling, cultivation, and design 
of the provenance trial. 


Sampling and traits 


To test for leaf morphological adaptation to 
climate, we identified a subset of 180 trees within 
the provenance test representing the extremes and 
median of multivariate maternal climatic conditions 
based on a principal component analysis of a 30-year 
(1971-2000) georeferenced baseline dataset of 23 
annual climate variables (Wang et al. 2016) at the 
source locations for the 672 maternal lineages. Six 
maternal lines were selected from among each of the 
highest, lowest, and median scores for the first 
principal component of all climate variables with 
the stipulation that no more than two maternal lines 
were selected from one locality (Fig. 1). The first 
principal component explained 43.8% of the variance 
in the climate data and was most strongly associated 
with variables involving temperature and seasonality 
(strongest positive factor loadings included mean 
annual temperature [0.28], degree-days > 5°C [0.28], 
and number of frost-free days [0.25]; strongest 
negative factor loadings included degree-days below 
18°C [-0.30], mm precipitation as snow [-0.27], and 
degree-days below 0°C [-0.25]). We haphazardly 
sampled five leaves from all progeny (n = 5) from 
each of the 18 focal maternal lineages within each of 
the two experimental sites (Chico and Placerville; 
because progeny mortality led to incomplete sam- 
pling, the total trees sampled = 141 and total leaves 
sampled = 707). We measured the area (mm”) and 
perimeter (mm) of each leaf from scanned images of 
fresh leaves using ImageJ (Schneider et al. 2012). 
Leaves were then dried to constant mass at 50°C, and 
dry biomass of each leaf was recorded to the nearest 
mg using an analytical balance. Measurements were 
used to calculate specific leaf area (area - dry 
biomass | in mm? - g') and leaf dissection (perim- 
eter - area’! in mm- mm”). 

To test the association between leaf morphology 
and photosynthesis, we selected trees to subsample 
within the Chico site using leaf trait data collected by 
MacDonald (2017), which included all progeny from 
all 672 maternal lines within both experimental sites. 
MacDonald’s (2017) study tested the effects of 
maternal climate on plasticity (or variation) in leaf 
traits, whereas the current study examined mean 
differences in traits among climate extremes as 
described by the PCA scores (see above). During 
2016, MacDonald (2017) measured morphological 
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traits on five leaves from each of the 3360 trees within 
each of the two common garden locations. Using this 
larger dataset, we selected 24 of the 672 maternal 
lines within Chico that had mean values of leaf 
dissection that were high (mean [SD] = 3.553 [0.157]), 
intermediate (2.628 [0.003]), or low (2.038 [0.056]) as 
reported by MacDonald (2017). Our objective with 
this approach was to sample trees representing a 
broad range of leaf dissection values across the 
experiment. By focusing on these lines, we assumed — 
and were able to test (see Results) — that variation in 
leaf morphology between years within a maternal line 
would be smaller than variation among lines. 

Using a portable infra-red gas analyzer (Li-Cor 
6400, LI-COR Biosciences, Lincoln, NE) we mea- 
sured carbon assimilation (A, pmol CO, -m~- s! 
on three leaves from each of a total of 28 progeny 
chosen among 17 of the 24 targeted maternal lines to 
maximize representation across expected leaf dissec- 
tion values, based on measurements from 2016 
(MacDonald 2017). Time constraints required for 
measurements with the instrument (e.g., internal 
equilibration required about 15 min for each leaf 
before measurements could be recorded) limited the 
number of trees we were able to sample. Measure- 
ments were collected across six days in 2018 (18, 19, 
26, 29 June and 9, 13 July) in the morning between 
0800 and 1200 h, prior to midday stomatal closure 
when photosynthetic rates fell to zero. The site was 
irrigated once during this period on 22 June 2018. 
Leaves were selected based on appearance (no 
apparent damage from herbivores). Internal chamber 
settings were as follows: photosynthetically active 
radiation (PAR) = 1500 pmol - s~'; stomatal ratio = 
1.0; flow = 500 pmol - s'; and reference chamber 
CO> concentration = 400 ppm. For leaves that were 
too small or deeply lobed to completely fill the area 
of the 2 cm? chamber, the area of leaf tissue inside the 
chamber was calculated manually by tracing the 
edges of the leaf inside the chamber with ink, 
removing the leaf from the tree, and later scanning 
an image of the leaf using a digital scanner. We then 
used ImageJ (Schneider et al. 2012) to estimate the 
measured area of each leaf, and gas exchange 
estimates were adjusted based on these values using 
the Li-Cor recompute utility. In addition, we used 
ImageJ to measure the area (mm) and perimeter 
(mm) of each leaf, which were then used to calculate 
leaf dissection as described above. 


Analysis 


Analyses were conducted using SAS version 9.3 
(SAS Institute, Cary, NC). Specific leaf area and leaf 
dissection were each compared across climate cate- 
gories (high, median, low) using mixed-model nested 
ANOVAs with the Kenward-Roger approximation 
for degrees of freedom. Climate category and local- 
ity(category) were included as fixed effects, and 
maternal family (category x locality) and individual 
tree (category x locality x family) were included as 
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TABLE 1. MIXED-MODEL ANALYSES OF VARIANCE OF SPECIFIC LEAF AREA (MM - G') AND LEAF DISSECTION (MM - 
MM ~) AMONG QUERCUS LOBATA LEAVES SAMPLED FROM 180 SAPLINGS (10 PROGENY FROM EACH OF 18 FAMILIES, 
ORIGINATING FROM 9 LOCALITIES) PLANTED IN TWO COMMON GARDENS. Sampled trees were selected to represent 
extremes and median of multivariate climatic conditions of maternal source trees (see Methods). F-tests of fixed effects and 
y° values from likelihood-ratio tests of random effects are shown. 


Test site Trait Effect type Source Fade gat OF 2 P-value 
Chico Specific leaf area Fixed Maternal climate category 0.615 65 0.5 
Locality (category) 1.14665 0.3 
Random Family (category x locality) 0.0 0.9 
Individual tree (category x locality x family) 163.9 <0.0001 
Leaf dissection Fixed Maternal climate category 0.025, 65 0.9 
Locality (category) 1.416. 65.4 0.2 
Random Family (category x locality) 0.0 0.9 
Individual tree (category x locality x family) Yom <0.0001 
Placerville Specific leaf area Fixed Maternal climate category 0.545 57.2 0.6 
Locality (category) 1.066, 57.3 0.4 
Random Family (category x locality) 0.0 0.9 
Individual tree (category x locality x family) 40.9 <0.0001 
Leaf dissection Fixed Maternal climate category 0.35 9.11 OT 
Locality (category) 1.666, 8.91 0.2 
Random Family (category x locality) 2.4 0.06 
Individual tree (category x locality x family) 43.2 <0.0001 


random effects. The significance of random effects 
was tested using likelihood-ratio tests (Littel et al. 
1996). To correct for heteroscedasticity, analyses 
were conducted using natural log-transformed values 
of leaf dissection, although means are reported below 
on the untransformed scale. Specific leaf area 
required no transformation to meet model assump- 
tions. Separate statistical models were used for data 
collected from each site (Chico and Placerville). 
Significant differences in leaf traits among climate 
categories were considered consistent with the 
hypothesis that leaf shape evolves in response to 
local climate. For the Chico dataset focused on 
carbon assimilation and leaf morphology, variation 
in A (carbon assimilation) was regressed against leaf 
dissection using a mixed-model nested ANCOVA. 
Sampling date and leaf dissection were included as 
fixed effects, and individual tree was included as a 
random effect. The sampling date x leaf dissection 
fixed effect interaction term did not significantly 
explain variation in A (P > 0.6), indicating that the 
effect of leaf dissection on A did not vary among 
sampling dates, so it was removed from the final 
model. Leaf dissection and A were both natural log- 
transformed prior to analysis to correct for hetero- 
scedasticity. A significantly positive relationship 
between leaf dissection and A was considered 
consistent with the hypothesis that greater leaf 
dissection is associated with higher photosynthetic 
rates. 


RESULTS 


Leaf dissection and specific leaf area varied widely 
among sampled leaves (coefficient of variation [CV] 
for specific leaf area = 28.4 and 19.0 for Placerville 
and Chico, respectively; for leaf dissection, CV = 


29.8 and 25.4 for Placerville and Chico, respectively). 
Nonetheless, neither trait varied significantly among 
categories of maternal multivariate climatic charac- 
teristics at either site (Table 1; Fig. 2). 

The subsampling of trees for measuring A was 
successful at capturing variation in leaf dissection 
within Chico. The CV of maternal mean leaf 
dissection values from the full dataset measured in 
2016 was 13.6 (MacDonald 2017), whereas the CV 
from our 2018 subsample of 18 maternal mean leaf 
dissection values was 19.6. In addition, maternal 
mean leaf dissection scores were positively correlated 
across the two datasets (rz = 0.71, P = 0.0009, n = 
18), validating our assumption that measurements of 
leaf dissection on trees in 2016 would predict 
measurements during 2018. Carbon assimilation 
increased significantly with leaf dissection among 
sampled trees (Table 2; Fig. 3; Slope = 0.86). In 
addition, estimates of A declined over the sampling 
period, possibly reflecting increasing air temperature 
and water limitation during the early summer 
sampling period (Table 2; Fig. 3). 


TABLE 2. MIXED MODEL ANALYSIS OF COVARIANCE 
TESTING THE EFFECTS OF LEAF DISSECTION (MM - MM ~) 
AND SAMPLING DAY ON CARBON ASSIMILATION (MOL 
CO, -M~-S') IN QUERCUS LOBATA GROWING IN A 
COMMON GARDEN EXPERIMENT AT CHICO. Twenty-eight 
saplings (17 families) were subsampled to maximize 
variation in leaf dissection (see Methods). F-tests of fixed 
effects and y° value from a likelihood-ratio test of the 
random effect are shown. 


Effect type Source Fat dat OF 7 P 

Fixed Leaf dissection Dil Sawa 0.03 
Sampling day 9.55 20.3 <0.0001 

Random Individual Tree 8.6 0.003 
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FIG. 2. Tukey’s box plots of specific leaf area and leaf dissection (perimeter - area) among leaves collected from Quercus 
lobata grown in common gardens in Chico (A and B) and Placerville, California (C and D) as a function of multivariate 
climatic characteristics of the maternal seed source (see Methods). See Table 1 for results of statistical tests. 
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Fic. 3. Carbon assimilation as a function of leaf dissection (perimeter - area | in mm - mm ~) in Quercus lobata sampled at 
the Chico common garden site. Sampling day was included in the model as a covariate. Carbon assimilation and dissection 
were log-transformed for analysis (see Methods). See Table 2 for results of statistical tests. 
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DISCUSSION 


Oak leaves are well known to harbor tremendous 
variation in shape. We explored two hypotheses 
associated with that variation. First, we tested the 
hypothesis that variation in leaf shape among 
populations is a consequence of adaptation to 
prevailing climatic conditions. Previous studies have 
reported that leaf morphological traits, including leaf 
dissection, vary with climate and other site charac- 
teristics (Ogaya and Penuelas 2007; Peppe et al. 2011; 
Royer 2012; Albarran-Lara et al. 2015; Riordan et al. 
2016; MacDonald 2017). If this hypothesis were true, 
we would expect variation in leaf morphology among 
trees in the common garden to be associated with the 
climatic conditions of the maternal source trees, 
because of inherited leaf morphological variations 
from the maternal trees. We instead found that the 
wide variation in leaf morphology observed was not 
associated with variation in climate, thus failing to 
support this hypothesis (Fig. 2). This may indicate 
that leaf dissection is principally influenced by the 
environment. Trait plasticity in response to climate 
may itself be adaptive in variable environments and 
could facilitate species persistence in response to 
climate change in circumstances where adaptive 
evolution is slow (Hoffmann and Sgro 2011; but 
see MacDonald 2017). On the other hand, pheno- 
typic variation in leaf morphology may be genetically 
controlled, but not associated with adaptation to 
climate. MacDonald (2017), in a much larger study 
of leaf morphology within the common garden, 
reported significant variation among family lines 
for leaf dissection and specific leaf area, which 
suggests that variation is genetically based. 

The extent to which leaf morphological variation 
in natural populations of Q. /obata is associated with 
climate is unclear, although the data herein may 
provide additional insight into this hypothesis. 
Albarran-Lara et al. (2015) examined phenotypic 
variation in Q. lobata leaves sampled from 127 adults 
distributed among three elevational transects and 
found substantial morphological variation, although 
climatic variation within transects was relatively 
small. Leaf thickness in the progeny of the sampled 
trees was positively correlated with maximum tem- 
perature of the site of maternal trees (Albarran-Lara 
et al. 2015). Other studies of oaks have reported that 
leaf morphology varies adaptively in response to 
climate (Abrams 1990; Cavender-Bares et al. 2004; 
Riordan et al. 2016; but see Royer et al. 2008). We 
suggest that study of phenotypic variation in field 
populations of valley oak across its range may 
provide additional insight into the validity of the 
hypothesis that leaf morphology is associated with 
climate. 

Second, we tested the hypothesis that variation in 
_ leaf shape has a functional consequence for plant 

physiology. We found that more dissected leaves had 
higher photosynthetic rates (Fig. 3), which is 
consistent with this hypothesis. Carbon assimilation 
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(A) measurements were collected during June and 
July, and A decreased over the period during which 
data were collected. California’s central valley, where 
Chico is located, has a very hot, dry season typically 
beginning in May; during the study no precipitation 
fell and temperatures routinely exceeded 35°C 
(University of California, Division of Agriculture 
and Natural Resources 2019). Even though the site 
was irrigated for one day during the time we collected 
A measurements (see Methods), the intensity of water 
stress experienced by the trees may have increased 
throughout this period. We avoided collecting A 
measurements for four days following the irrigation 
event to minimize its impact on our data. In any 
event, declines in A with increasing water stress are 
not uncommon (Chaves 1991; Osakabe et al. 2014), 
thus the temporal decline in A that we observed may 
have been a consequence of increased water stress 
experienced by trees. 

The positive relationship between leaf dissection 
and A may be a consequence of water relations. Prior 
research, for example, has found variation in leaf 
dissection in oaks to be negatively associated with 
hydraulic resistance (Siso et al. 2001). Similarly, 
Nicotra et al. (2008) found higher A and greater 
water loss in leaves of Pelargonium that were more 
dissected. In addition, other studies have suggested 
that a dissected morphology improves air exchange 
across leaf surfaces, leading to a more efficient 
transfer of heat (Gurevitch and Schuepp 1990). Thus, 
as early summer became increasingly hot, individual 
trees with more dissected leaves may have experi- 
enced reduced water stress, and consequently, these 
individuals may have been able to maintain higher 
photosynthetic rates under dry conditions than trees 
with less dissected leaves. Although we were not able 
to collect these data during the current study, this 
idea could be tested by concurrent collection of plant 
water potential and carbon assimilation measure- 
ments in trees with a range of leaf dissection values. 

This study provides evidence that variation in leaf 
shape found among valley oak individuals within the 
common garden is associated with photosynthetic 
rate. In addition, our results show that factors other 
than maternal climate influence variation in the leaf 
traits we measured. Some of these factors may be 
environmental, as the traits we measured are known 
to hold substantial phenotypic plasticity; MacDon- 
ald (2017), for example, found leaf dissection to be 
the most plastic of five leaf traits measured. Studies 
of other oak taxa also report plasticity in leaf 
morphology that is associated with functional 
responses to environmental variation (Wu et al. 
2016; Kusi and Karsai 2019). Within the Q. lobata 
common garden, genetically based variation in leaf 
morphology has also been reported (MacDonald 
2017). Oaks are known to harbor variation in leaf 
traits that involve adaptations to climate (Cavender- 
Bares 2004; Ramirez-Valiente et al. 2017a; Mead et 
al. 2019). Variation in leaf shape within valley oaks, 
whether plastic or genetic, appears to play a 
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functional role in mitigating environmental stress 
and climatic variation. Our results suggest future 
studies exploring the functional consequences of 
variation in leaf shape would be productive. 
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ABSTRACT 


We update and revise the serpentine affinity database of Safford et al. (2005), which documents levels of 
plant taxon association with ultramafic (“serpentine”) substrates in the California flora. The revised database 
reflects recent taxonomic changes consistent with the second edition of the Jepson Manual (TJM2) and more 
recent updates to the Jepson eFlora, and includes additional species that were not previously documented as 
serpentine associates. We also include serpentine-associated species that have been described since the 
publication of TJM2 and are not yet incorporated in the eFlora. A number of taxa were removed from the 
Safford et al. (2005) database due to new ecological information or taxonomic changes. As before, the new 
database presents information on rarity, geographic distribution, taxonomy, and lifeform. Based on our new 
database numbers and TJM2’s list of California endemic species, the percentage of California endemic full 
species that are +/— restricted to ultramafic substrates has risen to 14.7%. Of 255 total endemic taxa in our 
database (including strict and “broad” endemics), 148 (c. 60%) come from only ten plant families, 
concentrated mostly in one or two genera per family. The North Coast and Klamath Ranges continue to 
support more serpentine endemics than the rest of the State combined. The previous version of the database 
has been widely used in biodiversity research and conservation management, and we hope that this improved 


and updated version will prove as valuable. 


Key Words: California flora, edaphic, endemics, serpentine, ultramafic. 


Soils derived from ultramafic parent materials— 
often referred to as “serpentine’—support a great 
diversity of plant life in California. Serpentine soils 
are generally ecologically stressful environments 
because they are poor in plant macronutrients, have 
high magnesium : calcium ratios, and contain toxic 
heavy metals (Alexander et al. 2007). Serpentine 
environments support high levels of native plant 
diversity partly because they are less likely to support 
the (often non-native) dominant species that com- 
petitively exclude other plants on more fertile 
substrates (Brooks 1987). In addition, the stressful 
environment of serpentine leads to high levels of 
endemism. While serpentine soils occupy less than 
1.5% of California’s geographic area (Harrison et al. 
2000), a high number of California’s endemic plant 
species are partially or entirely restricted to serpen- 
tine soils (described below, also see Safford et al. 
2005; Safford 2011). Understanding the plant diver- 
sity associated with serpentine environments is 
important for facilitating research into basic ecolog- 
ical and evolutionary questions (of which there is a 
long history in California serpentine; Harrison and 
Rajakaruna 2011), as well as for informing manage- 
ment decisions involving landscapes, natural distur- 
bances, and conservation (Safford and Harrison 
2004, 2008). 

To this point, the serpentine endemism database 
of Safford et al. (2005) has provided the most 
comprehensive documentation of plant species affin- 


ities for serpentine substrates in California. This 
database assigns plant taxa a ranking representing 
their degree of specialization to serpentine soils based 
on herbarium records, botanical descriptions, field 
data, expert opinion, and published literature. The 
database has been widely used for a variety of 
purposes since its publication. The rankings have 
been cited by numerous ecological and evolutionary 
studies, such as research on the evolutionary origins 
of serpentine adaptation (Anacker et al. 2010), 
biogeographic considerations of plant diversity and 
endemism in the California flora (Harrison et al. 
2006; Burge et al. 2011), autoecological studies of 
serpentine species (Eckert 2006), and studies of 
mycorrhizae on serpentine soils (Southworth et al. 
2014). The serpentine affinity rankings have been 
integrated into the CalFlora (http://www.calflora. 
org) web interface—a database on California 
plants—where the rankings are displayed on the 
pages of serpentine-associated species. The serpen- 
tine database has also been used by land managers in 
conservation planning, environmental assessments 
and impact statements, and resource management 
plans. 


METHODS 


Our understanding of the systematics and ecology 
of the California flora has expanded substantially 
since the earlier version of this database was 
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TABLE 1. TAXA DROPPED FROM SERPENTINE AFFINITY DATABASE AFTER 2012. TJM2 (Baldwin et al. 2012) included 
many taxonomic changes and new taxa. These changes were incorporated into the database at the time of TJM2. This table 
includes taxa dropped from the list since that time. 


Species Family 
Chlorogalum purpureum var. reductum Agavaceae 
Hoover 
Cardamine pachystigma var. Brassicaceae 
pachystigma (S. Watson) Rollins 
Sedum obtusatum ssp. obtusatum A. Crassulaceae 
gray. 
Carex amplectens Steud. Cyperaceae 
Swertia fastigiata Pursh Gentianaceae 
Calochortus greenei S. Watson Liliaceae 
Ceanothus papillosus Torr. & A. Gray | Rhamnaceae 
var. roweanus McMinn 
Brodiaea pallida Hoover Themidaceae 
Brodiaea purdyi Eastw. Themidaceae 


Reason and source 


Taxon does not grow on serpentine soils (Kofron et al. 2013) 


Taxon is not a serpentine indicator; reanalysis of herbarium 
and collection records 

Taxon is not a serpentine indicator; reanalysis of herbarium 
and collection records 

Moved to Carex fracta Mack. by Baldwin et al. (2012), which 
is not a serpentine indicator 

Moved to Frasera umquaensis M. Peck & Applegate by 
Baldwin et. al (2012), which is not a serpentine indicator 

Taxon does not grow on serpentine soils; government agency 
(BLM, USFWS) reports on species ecology 

Taxon not recognized by TJM2 or Jepson eFlora, herbarium 
collections and CalFlora observations do not support taxon 
as serpentine indicator 

New collections show that only | of 5 population groups 
occurs on serpentine; R. Preston, UC Davis, pers. comm. 

Taxon subsumed within B. minor (Benth.) S. Watson, which 
is not a serpentine indicator (Preston 2006; Baldwin et al. 
2012) 


published. Here, we provide a revised version of the 
database (Appendix 1; Supplemental Table 1), with 
nomenclature updated to reflect The Jepson Manual, 
2nd edition (TJM2; Baldwin et al. 2012) and then 
again using the most recent Jepson eFlora (most 
recent version accessed November 18, 2019; Jepson 
Flora Project 2019). This contribution follows on 
Safford (2011), which provided a lay person’s 
summary of the status of the endemism database 
just as TJM2 was being prepared for press. The 
general patterns discussed in Safford (2011) are still 
valid, but it was important to incorporate the 
published results from Baldwin et al. (2012) and the 
most recent Jepson eFlora (Jepson Flora Project 
2019), and to provide a completely revised version of 
the database itself, which was not part of Safford 
Cory. 

In the current database (Appendix 1), we have 
added numerous additional species that were not 
included on the original 2005 list, many of which are 
recently described serpentine-associated species. We 
also removed seven taxa from the list that we 
constructed after TJM2 (Table 1) because recent 
taxonomic changes and/or further study has indicat- 
ed they are not strongly associated with serpentine 
substrates. To revise the original list, we first updated 
nomenclature to current Jepson eFlora names (as of 
November 18, 2019; Jepson Flora Project 2019). We 
then solicited feedback on the list from numerous 
expert California botanists and requested suggestions 
for additional species that should be included. We 
also conducted a literature review for recently 
described species that have affinities for serpentine. 
We then conducted further research on all species 
that were identified by experts or through the 


literature review as possibilities for addition to or 
removal from the list, as described below. 

We have ranked new additions to the list using the 
same semi-quantitative scale as the original database 
(see Safford et al. 2005 for details), with the 
difference that we have reduced all rankings >6 to 
6, as those resulted from the addition of CalFlora 
information (which at the time could be considered 
independent of our other sources). Today CalFlora 
directly references the serpentine affinity database in 
its species descriptions. Rankings of new additions to 
the list presented here are based primarily on 
herbarium records (mostly accessed via the Consor- 
tium of California Herbaria (CCH1 2019; CCH2 
2019), habitat descriptions by Baldwin et al. (2012) 
and the Jepson eFlora (Jepson Flora Project 2019), 
CalFlora observation records (CalFlora 2019), the 
primary literature, and in some cases, direct input 
from botanical experts. In the CCH herbarium 
records, we accessed at least 150 collection records 
for each taxon (except where records totaled <150, in 
which case we accessed all records), using random 
selection where necessary to ensure no bias in 
regional coverage. Where CCH or CalFlora records 
were <5, we combined the two and treated them as a 
single information source. We used the Berkeley 
Mapper (linked to the CCH website), Google Earth, 
our field experience, and various geologic and 
jurisdictional (e.g., US Forest Service) maps to 
visualize the locations of herbarium records and to 
determine, when possible, whether botanical records 
lacking geological information were found on or off 
of ultramafic substrates. Taxa with >50% of usable 
records (“usable”’ = the records identified the 
substrate, or we were able to determine the substrate 
by another fashion) from serpentine made their way 
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TABLE 2. TAXA INCLUDED IN THE SERPENTINE AFFINITY LIST THAT ARE NOT IN THE CURRENT JEPSON EFLORA (AS OF 
NOVEMBER 18, 2019). These include recently described taxa as well as taxa that are widely recognized by botanists in 
California, but are not included in the current eFlora because of taxonomic uncertainty. 


Taxon 


Brodiaea rosea ssp. rosea (Greene) Baker 
Calystegia vanzuukiae Brummitt & Namoff 


Cardamine pachystigma vat. dissectifolia (S. Watson) Rollins 


Claytonia serpenticola T.R. Stoughton 
Dicentra formosa ssp. oregana (Haw.) Walp. 


Erythranthe willisii G.L. Nesom 

Hesperolinon sharsmithiae R. O'Donnell 

Sedum eastwoodiae (Britton) A. Berger 

Sedum flavidum B.L. Wilson & Zika 

Sedum kiersteadiae B.L. Wilson & R.E. Brainerd 
Sedum patens Zika 

Sedum rubiginosum Zika & B.L. Wilson 
Vaccinium coccineum Piper 


onto the final list. Rankings of new list additions 
(e.g., recently described species) should be used 
cautiously since many of them are based on limited 
data, and we suggest that all species ranked based on 
three or fewer sources are in need of further study. 
All species in the updated list follow the nomencla- 
ture of the current (November 2019) Jepson eFlora 
except for 13 taxa (Table 2), most of which are 
recently described. 

We emphasize that we did not carry out a 
complete search of new collection records for every 
taxon since 2005. This was only carried out for taxa 
where recent publications, taxonomic changes, or 
personal communications from California botanists 
led us to believe that such a search might change the 
status of the species. For new additions to the list, we 
no longer required three separate sources of infor- 
mation supporting their affinity for serpentine if 
evidence for their affinity to serpentine appeared 
strong (for example, in the herbarium records). 

Finally, it is important to note that for species for 
which all described infraspecific taxa (subspecies, 
varieties) exhibit the same level of serpentine affinity, 
we only include the full species name in our database. 


RESULTS AND DISCUSSION 


As in Safford et al. (2005) and Safford (2011), we 
define serpentine endemic taxa as those that scored at 
least 4.5 in our ranking process. Put another way, our 
class of “serpentine endemics” includes strict endem- 
ics (more than c. 95% of recorded occurrences on 
ultramafics) and broad endemics (c. 85-94%). The 
updated serpentine affinity database (Appendix 1) 
contains 255 endemic taxa (strict plus broad endem- 
ics), 193 of which are full species (this includes five 
species for which all of the subspecies are listed in the 
database due to differing levels of serpentine affinity: 
Arctostaphylos bakeri Eastw., A. montana Eastw., 
Calystegia collina (Greene) Brummitt, Cirsium fonti- 


Source 


Preston 2013 

Brummitt and Namoff 2013 

R. Preston and L. Janeway, pers comm and CCH 
and CalFlora records 

Stoughton 2017 

Baldwin et al. 2012; L. Hoover, U.S. Forest 
Service, pers. comm. 

Nesom 2017 

O’Donnell 2006 

Zika et al. 2018 

Zika et al. 2018 

Wilson et al. 2014 

Zika et al. 2018 

Zika et al. 2018 

L. Hoover and J. Nelson, U.S. Forest Service, 
pers comm. 


nale (Greene) Jeps., and Lessingia micradenia 
Greene). Thus, serpentine endemic species account 
for about 14.7% — more than 1/7 — of the endemic 
species in California (based on 1315 endemic species 
in Baldwin 2012). This is a notable increase from the 
12.5% reported by Safford et al. (2005), and is due to 
additions of newly described species and a decrease in 
the number of California endemic species according 
to Baldwin et al. (2012). 

Taxonomic revisions have led to some changes 
from Safford et al. (2005) in the ordering of family 
and genera with the most serpentine endemic taxa. 
Asteraceae remains the family with the most serpen- 
tine endemics, with 25 strict endemics, 13 broad 
endemics, and 106 total serpentine-associated taxa 
(Table 3). Liliaceae, which historically contained the 
second most serpentine endemics of any family 
(considering both strict and broad endemics), has 
been split into several smaller families, many of 
which are further down the current list, including 
Alliaceae — also in the top ten families, Agavaceae, 
Themidaceae, and so on. As a result, Brassicaceae is 
now the second most endemic-rich family (Table 3). 
Scrophulariaceae was also split in Baldwin et al. 
(2012) into numerous plant families, and it has 
dropped from a fifth place ranking in the Safford et 
al. (2005) rankings to a complete absence from the 
list now. Families derived from Scrophulariaceae are 
still important components of serpentine endemism 
in California however: for example, Orobanchaceae 
currently occupies sixth place on the list, and 
Phrymaceae and Plantaginaceae are further down 
the list. Eriogonum is now the most endemic-rich 
genus, followed by Streptanthus (Table 4). 

As discussed by Safford et al. (2005; see their Figs 
1 and 2), the distribution of taxa among plant 
families and genera is very different when comparing 
serpentine endemics and California endemics. Fam- 
ilies which support many more serpentine endemic 
taxa than expected include Linaceae, Polygonaceae, 
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TABLE 3. NUMBERS OF SERPENTINE ENDEMIC AND NEAR 
ENDEMIC TAXA BY FAMILY. | Strict endemics. * Strict 
endemics plus broad endemics (“endemic taxa” in this 
paper). * Strict and broad endemics plus ‘“‘near endemic” 
taxa (taxa transitional from strong indicators to broad 
endemics). * Total taxa in our database, including weak 
indicators and taxa transitional from indifferent to weak 
indicators. Families ranked by numbers of endemics (strict 
+ broad). 


Serpentine affinity score 


Family 55" >4.5° >35?. <Total* 
Asteraceae 25 38 46 106 
Brassicaceae 19 24 28 38 
Polygonaceae IFIR 19 yp 4] 
Liliaceae 5 iB 16 38 
Polemoniaceae 9 10 1] 21 
Orobanchaceae yi 10 10 19 
Apiaceae 7 9 12 33 
Crassulaceae 8 9 9 16 
Linaceae 8 9 9 14 
Ericaceae 5 8 10 15 
Fabaceae 4 a: i 25 
Alliaceae 4 7 11 24 
Lamiaceae 6 7 9 1 
Boraginaceae 6 6 1] 24 
Caryophyllaceae 5 6 8 18 
Campanulaceae 3 5 8 13 
Rhamnaceae 4 >) i 14 
Onagraceae 3 5 7 13 
Cyperaceae 3 5 y 12 
Convolvulaceae l 5 6 7 
Phrymaceae 4 4 6 8 
Rubiaceae 3) 4 4 8 
Montiaceae 0 3 6 17 
Rosaceae ] 3 6 1] 
Plantaginaceae ] 3 p LS 
Poaceae ] 3 3 18 
Iridaceae 2 3 3 5 
Violaceae 0 2 3 7 
Cupressaceae 0 2 3 6 
Malvaceae yp 2 js 7 
Agavaceae 2 2 2 5 
Salicaceae 2 2 2 3 
Gentianaceae p 2 2 2 
Garryaceae l 2 2 2 
Themidaceae 0 1 5 8 
Ranunculaceae l f 3 ii 
Papaveraceae 0 l I s) 
Berberidaceae 1 1 l 4 
Pteridaceae 0 1 1 4 
Fagaceae | ] ] é) 
Orchidaceae 0 1 1 3 
Dryopteridaceae l ] I 2 
Lentibulariaceae ] l 1 1 
Apocynaceae ] l ] | 
Verbenaceae 0 1 1 l 
Melanthiaceae 0 0 2 4 
Pinaceae 0 0 1 6 
Parnassiaceae 0 0 l v4 
Saxifragaceae 0 0 | ] 
Sarraceniaceae 0 0 l 1 
Tecophilaeaceae 0 0 0 ] 
Primulaceae 0 0 0 l 
Polygalaceae 0 0 0 ] 
Cistaceae 0 0 0 l 
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TABLE 4. GENERA WITH MORE THAN THREE TAXA 
ENDEMIC TO SERPENTINE (MEAN AFFINITY SCORE > 4.5). 


Genus Family Number of taxa 
Eriogonum Polygonaceae 19 
Streptanthus Brassicaceae 18 
Allium Alliaceae 1] 
Arctostaphylos Ericaceae 9 
Hesperolinon Linaceae go 
Carex Cyperaceae 7 
Lomatium Apiaceae of 
Erythranthe Phrymaceae 6 
Navarretia Polemoniaceae 6 
Packera Asteraceae 6 
Phacelia Boraginaceae 6 
Sedum Crassulaceae 6 
Calystegia Convolvulaceae 5 
Campanula Campanulaceae 5 
Cordylanthus Orobanchaceae 5 
Erigeron Asteraceae 5 
Fritillaria Liliaceae 5 
Monardella Lamiaceae 5 
Boechera Brassicaceae 4 
Calochortus Liliaceae 4 
Castilleja Orobanchaceae + 
Ceanothus Rhamnaceae 4 
Cirsium Asteraceae 4 
Galium Rubiaceae 4 
Harmonia Asteraceae 4 
Horkelia Rosaceae + 
Lessingia Asteraceae 4 


Brassicaceae, Apiaceae, Caryophyllaceae, and Lilia- 
ceae. Genera that follow the same pattern of notable 
enrichment in serpentine endemics include Streptan- 
thus, Minuartia, Cordylanthus, Packera, and Hesper- 
olinon. 

Of the five geographic regions identified by 
Safford et al. (2005), the North Coast region remains 
the geographic leader in number of serpentine 
endemic taxa (127), and the Klamath region supports 
the most serpentine endemics that are restricted to a 
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FIG. 1. Geographic distribution of serpentine endemic 


taxa (mean score > 4.5) in California. Dark areas of bars 
represent taxa restricted to a single geographic region, and 
light areas represents taxa that also occur in other regions. 
Geographic regions are from Baldwin et al. (2012). 
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FIG. 2. Lifeforms of serpentine endemic and near endemic 
taxa (mean score >= 3.5) in California. 


single region (60; Fig. 1). Together, the North Coast 
and Klamath regions support 146 serpentine endemic 
taxa that are found nowhere else in the state. See 
Safford (2011) for a figure and short discussion 
regarding the elevational distribution of serpentine 
endemic taxa in California. In short, most taxa grow 
in habitats below 1000 m elevation and few are found 
above 1500 m, primarily due to the relative rarity of 
ultramafic outcrops above that elevation. 

Among the endemic taxa, 216 are dicots, 35 are 
monocots, two are gymnosperms, and two are 
pteridophytes (Appendix 1). The distribution of life 
forms is shown in Fig. 2. Serpentine endemic taxa in 
California are primarily perennial forbs (153 of 255 
taxa, including three that can also be annuals), most 
of the rest are annual forbs (74 taxa, including the 
same three taxa); less than 10% of the endemic taxa 
in our database are woody species (Fig. 2). 

192 (77%) of the endemic taxa in our database are 
listed in the CNPS Inventory of Rare and Endan- 
gered Plants (CNPS 2019). 108 serpentine endemic 
taxa are assigned California Rare Plant Rank 1B, 
which includes Federal and State threatened and 
endangered taxa as well as those taxa that the 
California Native Plant Society believes warrant 
listing (Appendix 1). 

The richness of serpentine endemic taxa in 
California is astonishing, given the limited area of 
the state that supports ultramafic rocks and soils. 
Yet, amazingly, new species continue to be described 
from ultramafic soils in California. Only 15 years 
after embarking on this exercise for the first time, we 
expanded the list of California serpentine endemics 
by 3.7%. Given the outsized contribution made to 
California’s natural riches by ultramafic substrates 
and their biota, we hope and expect that botanical 
and geoecological exploration will continue to 
expand the California serpentine flora, as well as 
our understanding of the factors that make it so 
diverse. 
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ABSTRACT 


Successful restoration of California terrestrial ecosystems requires knowledge of seed germination rates for 
a diversity of native species. We quantified greenhouse germination of 25 native perennial forbs, shrubs, and 
grasses collected coastally near Santa Cruz, CA for restoration at the Younger Lagoon Reserve and nine 
annual forb species from interior grasslands at Pinnacles National Park. Seeds of most species were collected 
in multiple years. Twenty-seven of our 34 study species had at least 25% germination without complicated 
germination triggers. All but two of the 13 species that were collected in >3 years and germinated showed 
significant interannual variation in germination; seven species had a trend of decreasing percent germination 
with increasing age. These species- and collection-year specific results indicate a necessity to test seed stock 
when possible prior to seeding to increase restoration success, and highlight the importance of following best 


practices for seed collection and storage. 


Key Words: central California, coastal prairie, grassland, restoration, sage scrub, seed germination, seed 


longevity. 


Sage scrub and grassland ecosystems are critically 
threatened in California (Hoekstra et al. 2005) due to 
land conversion, non-native invasive species, and 
altered disturbance regimes (Vasey and Holl 2007). 
The long-term recovery of California native ecosys- 
tems depends on developing best management 
strategies for their restoration, in particular propa- 
gating and planting native seed to ensure conserva- 
tion of biodiversity and genetic diversity at the 
project site and increasing site resistance to invasive 
species (Stromberg et al. 2007; Vasey and Holl 2007). 
Restoration requires species-specific information on 
seed germination rates and seed longevity data for 
the huge variety of California native species. 

Seed germination rates are determined by many 
variables including seed age, storage practices, 
maternal effects, climate, and site conditions (Rajjou 
et al. 2008; Baskin and Baskin 2014). Past studies of 
orthodox seeds, i.e., seeds that survive drying and/or 
freezing events in storage, have shown that many 
species retain viability for decades, whereas seed 
viability decreases rapidly over time for other species 
(Priestley et al. 1985; Walters et al. 2005; Fenner 
2017). Furthermore, a number of past studies have 
found reduced seed longevity over time when seeds 
are exposed to conditions with higher moisture and 
temperature (Dickie et al. 1988; Brown and Briggs 


1991; Nguyen et al. 2012). Maternally expressed 
genetic factors strongly influence seed development 
and resulting seed mass, which both influence 
resulting germination (Wulff 1995; Wang et al. 
2016). Likewise, site conditions (e.g., soil moisture 
and nutrients) can affect seed germination (Fenner 
2017). 

This study aimed to quantify the percent germi- 
nation of 34 native California grassland and coastal 
species that are commonly used in restoration 
projects in central California, most of which were 
collected in multiple years. Our results supplement 
existing information on seed germination require- 
ments and rates (e.g., Emery 1988; Gulmon 1992; 
Keeley and Fotheringham 1998; Wall and Macdon- 
ald 2009) to better enable restorationists to collect 
and store high quality seed, and in turn reduce 
project costs and increase plant establishment (Bar- 
ton et al. 2016). 


MATERIALS AND METHODS 


This study was conducted with seed collected for 
restoration at Younger Lagoon Reserve (YLR) and 
Pinnacles National Park. YLR is a 29-hectare area of 
protected land (36.95370°N, —122.06654°W) adjacent 
to the Pacific Ocean that belongs to the University of 
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California Natural Reserve System and is dominated 
by coastal scrub and grassland habitat. Much of the 
site is currently being restored following a long 
history of agricultural land use. Seeds for the YLR 
study site were collected in Santa Cruz County along 
State Route 1 at Afio Nuevo State Park (37.11931°N, 
—122.30765°W), Wilder Ranch State Park 
(36.96080°N, —122.08345°W), Coast Dairies State 
Park (36.98439°N, —122.15570°W), Scaroni Farms 
(36.97800°N, -—122.13815°W ), 4 mile Beach 
(36.96633°N , -122.12272°W), and the University of 
California at Santa Cruz (UCSC) at West Marshall 
Field (37.01869°N, —122.07562°W). 

Pinnacles National Park (PNP) is a federally 
designated U.S. National Park about 65 km inland 
from the coast which includes chaparral, oak 
woodlands, and grassland ecosystems. Park staff 
are restoring a nearly flat 80.9-hectare degraded 
grassland near the east entrance to the park 
(36.49066°N, 121.18252°W). Seeds were collected at 
many sites within the boundaries of PNP. 

We tested germination rates for 25 native 
perennial grass, forb, and shrub species that are 
commonly used for restoration at YLR and 9 native 
annual forb species that are being reintroduced to 
PNP (Appendix 1). Sixteen species from YLR were 
collected in 3—7 consecutive years (Appendix 1). 
Nine species from YLR and nine species from PNP 
were collected in 1 or 2 years (Appendix 1). 
Nomenclature throughout follows Jepson Flora 
Project.(2019), 

Seeds were collected from a minimum of 50 plants 
and at multiple sites, when possible, to ensure 
genetic diversity, following standard seed collection 
recommendations (Maschinski et al. 2012). The 
relative amount of seed of individual species 
collected from different sites varied annually de- 
pending on seed abundance at a given site. Seed 
collections from a given year were usually combined 
into one collection container after cleaning to 
conserve storage space. All seeds were initially 
stored in paper bags and air dried for 6—12 months. 
Seeds were then cleaned manually by passing seeds 
through a minimum of two mesh screens (0.074—3.8 
mm) until the target seed was completely free of 
chaff. After cleaning, YLR seeds were primarily 
stored at ambient temperature, except for three 
species that were stored at —2°C (Appendix 1); all 
collections of individual species were stored at the 
same temperature. YLR species were stored in one- 
gallon plastic bags or paper envelopes (depending 
on volume), except for a few very-small seeded 
species in glass jars; most collections of the same 
species were stored in the same container type. All 
seeds from PNP were stored at ambient temperature 
in paper envelopes. 

We conducted three rounds of germination testing 
at the Jean Langenheim Greenhouses, University of 
California, Santa Cruz. All PNP seed testing started 
on 13 January 2017; a small set of YLR seed was 
tested on 3 February 2017; and most YLR seeds, 
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including species that were collected in >3 years, 
were tested starting 11-13 January 2018. All collec- 
tions of the same species were tested in the same 
round, and each species was only tested during one 
round. For each species X collection year combina- 
tion, 100 seeds were divided equally between four 
10.16 X 10.16 X 15.24-cm pots filled with Pro-mix HP 
mycorrhizae soil. In each pot, 25 pre-cleaned, pre- 
counted seeds were placed into a 5 X 5 seed grid to 
avoid aggregation. We placed a light covering of soil 
over the seeds, approximately 1.5X the seed length. 
Pots were maintained at ambient temperature in a 
greenhouse receiving some indirect and some direct 
light. Pots were misted twice daily by the greenhouse 
irrigation system and watered additionally when the 
soil surface was dry to the touch. We monitored 
germination weekly for 11 weeks and removed all 
visible germinants from the soil with forceps cleaned 
with 70% alcohol solution to reduce microbial 
contamination. 

For 2017 PNP seed testing, Amsinckia menziesii 
(Lehm.) A.Nelson & J.F.Macbr., Calandrinia men- 
ziesii (Hook.) Torr & A.Gray, Castilleja exserta 
subsp. venusta (A.Heller) T.I.Chuang & Heckard, 
Clarkia unguiculata Lindl., Lasthenia gracilis (DC.) 
Greene, and Madia sativa Molina were cold stratified 
at 2.78°C for 2—3 weeks (Emery 1988). For 2017 YLR 
seed testing, all seeds were initially placed in cold 
stratification without light at 2.78°C for five weeks 
and then moved to ambient temperature. No species 
were stratified in 2018 except for Sisyrinchium bellum 
S.Watson, which was stratified at 2.78°C for 6 weeks 
pre-sowing (Emery 1988). 

ANOVA was used to evaluate whether the percent 
germination varied significantly among collection 
years for species that were collected in >3 years. 
When the effect of collection year was significant (a = 
0.05), we conducted Tukey’s multiple comparison 


tests to compare individual years. 


RESULTS 


Germination rates varied greatly across species 
and collection year with 27 of the 34 species (76%), 
including all the annual species from PNP, exhibiting 
at least 25% germination in one or more collection 
years. No additional treatments were used to 
stimulate germination aside from stratification of 
several species as mentioned previously (Fig. 1, Table 
1). Artemisia californica Less., Frangula californica 
(Eschsch.) A.Gray, and Symphyotrichum chilense 
(Nees) G.L.Nesom did not germinate, and Carex 
gynodynama Olney, Danthonia californica Bol., and 
Sisyrinchium bellum had <3% germination in all 
collection years. Of the 13 species that were collected 
over a period of >3 years that germinated, all but 
two showed significant interannual variation, and 
seven (54%) showed a trend of decreasing germina- 
tion with increasing year since collection (Fig. 1, 
Appendix 2). 
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TABLE 1. PERCENT SEED GERMINATION OF ALL SPECIES COLLECTED IN 1 OR 2 YEARS. Values are mean germination 


percentages + 1SD for n= 4. 


Species Test year 
Amsinckia menziesii 2017 
Calandrinia menziesii 2017 
Castilleja exserta ssp. venusta 2017 
Chaenactis glabriuscula 2017 
Clarkia unguiculata 2017 
Delphinium sp. 2017 
Lasthenia gracilis 2017 
Lupinus bicolor 2017 
Madia sativa 2017 
Cirsium quercetorum 2018 
Danthonia californica 2017 
Erigeron glaucus 2017 
Juncus mexicanus 2018 
Lupinus variicolor 2018 
Scrophularia californica 2017 
Stipa pulchra 2017 


DISCUSSION 


Our results show that the vast majority of species 
studied had >25% germination without requiring 
seed pretreatment to break dormancy, despite the 
fact that many native species in California have 
complex germination triggers to ensure that germi- 
nation only occurs when microclimatic conditions 
are favorable, such as sufficient rainfall and light 
(Keeley 1987; Emery 1988; Vaughn and Young 
2010). While it is recommended to test viability of 
non-germinating seeds through a cut or crush test or 
tetrazolium chloride (Baskin and Baskin 2014), this 
was not possible for several of the tiny-seeded species 
we studied. A lack of appropriate germination 
triggers was likely the cause of little to no germina- 
tion observed for Sisyrinchium bellum and Frangula 
californica in all collection years. Emery (1988) 
suggests 1.5 months stratification for newer seeds 
and no treatment for 3—6 year old seeds for S. bellum, 
and reports that stored F. californica seeds require at 
least three months of stratification. Also, some 
species may require two or more years of after- 
ripening to allow for gas exchange, water penetra- 
tion, and neutralization of inhibitory chemicals to 
initiate germination events (Emery 1988; Baskin and 
Baskin 2014), a pattern consistent with the germina- 
tion of Eriophyllum staechadifolium Lag. Likewise, 
Holl (unpublished data) observed that S. bellum 
germinated in the second year after sowing in 
another grassland restoration study near YLR. 


Site Seed age (yr) Germination (%) 
PNP l 84.0 + 5.7 
5 76.0 + 15.0 
PNP 1 43.0, + DAG 
PNP 5 43.0 + 6.0 
PNP 1 79.0 + 14.4 
PNP | 60.0 + 19.9 
5 83.0 + 11.9 
PNP | 94.0 + 5.2 
PNP l 37.0 + 6.0 
5 26.0 + 6.9 
PNP | 91.0 + 3.8 
5 O70: = 83 
PNP | 85.0 + 3.8 
5 58.0 + 14.8 
YLR l 53.0 = Fn 
2 305238 
YER l 1.0 + 2.0 
YLR 2 7.0 + 3.8 
3 18.0 + 8.3 
YLR 2 22eG = OS 
2 34.0 + 9.5 
TLR 3 SPO. = 16:0 
5 50.0 2 
YER 2 57.0 = 6.0 
5 74,03 3072 
YLR 2 23.0 + 8.2 
6 47.0 + 11.9 


Percent germination varied strongly by collection 
year. One explanatory factor could be seed age, 
consistent with previous studies, as approximately 
half of the species that were collected in three or 
more years showed a decrease in percent germination 
over time (Gulmon 1992; Tielb6rger and Petru 2010; 
Barton et al. 2016). Our results showing that Bromus 
sitchensis var. carinatus (Hook. & Arn.) R.E.Brai- 
nerd & Otting. and Elymus glaucus Buckley (both 


Poaceae) seeds are short lived when stored at room 


temperature are consistent with (Dremmen 2003) 
who reported seeds of those species only last around 
three years. Three of the four Asteraceae species we 
tested decreased in germination over time, and 
Artemisia californica seed that was three to five years 
old had low germination in all years. Observations 
from restoration efforts at the University of Cal- 
ifornia, Santa Barbara concur that stored Asteraceae 
seeds tend to be short-lived (J. Luong, University of 
California, Santa Cruz, personal communication). 
A number of other factors besides seed age could 
also explain the effect of collection year on germina- 
tion, particularly given that the percent germination 
of four species (Elymus glaucus, Hordeum brachyan- 
therum Nevski, Horkelia californica Cham. & 
Schltdl., and Prunella vulgaris L.) varied by collection 
year but showed no obvious directional trend. Past 
studies have shown that the rate of seed germination 
depends on initial seed quality and health of the 
parent population upon collection, which inevitably 
varies annually (Rajjou et al. 2008; Kochanek et al. 
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2011; Nguyen et al. 2012), often based on rainfall in 
the prior year (Kochanek et al. 2010; Zani and 
Miller 2017). 

Collection practices may also partially explain the 
mixed germination trends we observed. Variation in 
local site conditions of the mother plant has been 
shown to explain germination variation in past 
studies (Gulmon 1992; Tielb6rger and Petru 2010). 
However, since PNP and YLR seed collections from 
different locations were pooled into one collection 
container to conserve space and labeled only with the 
collection year, we cannot tease out the separate 
effects of seed age and site on germination. Further- 
more, the majority of seeds used in this study were 
stored in ambient temperature and moisture condi- 
tions, whereas best management practices call for 
storing seed at —20°C and low moisture conditions to 
maximize seed longevity (Linington and Manger 
2014). 


RECOMMENDATIONS 


Our results suggest several recommendations for 
collecting, storing, and testing restoration seed 
stocks. First, in light of the species-specific rates of 
seed germination and the strong effect of collection 
year on germination, we recommend that restora- 
tionists test seed stock periodically to calibrate 
appropriate seeding rates and increase restoration 
success. Second, as a general rule of thumb, seeds 
should be used for restoration within a couple of 
years of collection to avoid potential reduced 
germination due to seed aging (Brown and Briggs 
1991). Third, our results highlight the importance of 
collecting seed in multiple years and at multiple 
locations to ensure the greatest genetic diversity 
possible and to increase germination success (Emery 
1988; Brown and Briggs 1991; Maschinski et al. 
2012). Fourth, best management practices should be 
followed, such as labeling collections from different 
locations separately and cold storing seeds that will 
be held for multiple years (Maschinski et al. 2012). 
To date, most of the knowledge of seed germination 
triggers and longevity for California native species 
amassed by horticulturalists and restoration practi- 
tioners have not been widely disseminated. Hence, 
our final recommendation is for enhanced sharing of 
written documentation of successful and failed 
strategies in seed collection, storage, and sowing 
practices, which will enhance the success of native 
plant restoration throughout California. 
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Book Review 


Field Guide to the Grasses of 
Oregon and Washington. By 
CINDY TALBOTT ROCHE, 
RICHARD E. BRAINERD, 
BARBARA L. WILSON, 
NICK OTTING, AND ROB- 
ERT C. KORFHAGE. 2019. 
Oregon State University 
Press, Corvallis, Oregon. 
460 pp. ISBN 978- 
0870719592 (paperback). 
Price: $35.00, paperback 
edition. 


This field guide to all of the grasses (except 
bamboos) that grow outside cultivation in Oregon 
and Washington includes 376 taxa. The book has five 
sections: an introduction (14 pages), keys to genera 
and species (75 pages), and species accounts with 
color photos, description and distribution map (346 
pages), followed by a glossary and index. Why am I 
reviewing it for California botanists? Because Cal- 
ifornia botanists will find this book a valuable 
supplement to our own field guides to grasses. Many 
of the grasses that grow in Oregon and Washington 
also grow in California, ranging from widespread 
native and introduced species to even some rare 
species. You would not expect to find California’s 
rare endemic species, Agrostis blasdalei Hitche., in 
the field guide, but much to our surprise, it was 
recently found on the southern Oregon coast. (Of 
course, it is no longer a California endemic). Some 
California serpentine and alpine and subalpine 
grasses are also well described and illustrated in this 
guide. It does not include all California grasses; 
California has habitats not well represented in the 
Pacific Northwest, such as vernal pools, chaparral, 
and hot southern deserts, and, therefore, supports a 
greater diversity of species. Nevertheless, no treat- 
ment of California grasses is better illustrated or as 
up-to-date as this field guide (other than the Flora of 
North America; Barkworth et al. 2003, 2007). 

People have developed many ways to use field 
guides to identify a grass in hand: some go right to a 
key, while others leaf through pictures looking for 
something similar. Most of us prefer to ask someone 
who knows grasses, but let’s say we’re on our own 
with only a field guide. Where to start, what features 
to examine first, how to tell that it is indeed a grass 
and not a rush or sedge? The introduction meets our 
needs at this initial stage. The authors, whose 
extraordinary talents brought this exceptional field 
guide to fruition, note that many plants with 
common names that include the word “grass” are 


not grasses; for example, bear grass (Yerophyllum 
tenax (Pursh) Nutt.) is not a grass. Then, they 
distinguish grasses from other graminoids, like 
sedges and rushes, describing the differences in stem 
and flower characters. For example, sedge flowers are 
protected by one bract while grass flowers are 
enclosed by two bracts: the lemma and palea. The 
introduction also includes a section entitled “Grass 
Structure and Vocabulary,” which is as important as 
the emphatic “READ ME” next to the title indicates. 
You can’t use the key if you don’t know the 
terminology used in grass anatomy. You'll return 
frequently to these pages as you work the key and in 
time you will become fluent in the vocabulary. The 
vegetative structures are well described, often with 
practical suggestions. For example, flowers (also 
known as florets) are arranged in spikelets that are 
subtended by two additional bracts, called glumes. 
Spikelet compression can be detected by laying a 
spikelet on a tabletop. If both glumes are visible, it is 
laterally compressed. If only one glume is visible, it is 
dorsiventrally compressed. If it takes a variety of 
positions, it 1s probably not compressed; it is 
cylindrical or terete. In my opinion, the authors 
could have beefed up this section on structure and 
vocabulary with more drawings, or they could have 
illustrated the glossary. Since this section of the guide 
should address the needs of the beginner, it would 
have been useful to include definitions and, prefer- 
ably, illustrations of a// the technical words used in 
the key and species treatments. I selected a couple of 
species treatments randomly to test this notion and, 
although most of the descriptive words are defined 
either in the “Vegetative Structures” section or in the 


glossary, the very first species treatment I turned to, 


on page 148, used a word undefined anywhere: 
“scurfy.” To address this deficiency, I will paste an 
illustrated glossary from another source on the inside 
back cover of my copy. 

Well-written keys, like those here, will deliver the 
answer 90% of the time, but they aren’t perfect nor 
are their users. Characters in keys usually describe 
mature, complete plants, but specimens are frequent- 
ly incomplete, immature, or senescent; the authors 
thus recommend collecting multiple specimens. If 
you reach a dead end following one lead, go back and 
follow the other lead. If that doesn’t work, go back 
one step more and try again. The authors of the keys 
are veterans of the Carex Working Group, who 
published a field guide to sedges for the Pacific 
Northwest and also wrote the graminoid keys for the 
Flora of Oregon and the Flora of the Pacific 
Northwest (Meyers et al. 2015; Hitchcock and 
Cronquist 2018). They note that users will learn 
from experience that some characters are generally 
more reliable than others, that some grasses favor 
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particular habitats, and how to visualize the charac- 
ters in the descriptions. The emphasis on learning 
over time suggests that keying grasses requires 
patience, that even experienced naturalists find 
grasses challenging. 

The descriptions of species are generous in 
technical detail. Photographs of a whole plant, its 
habitat, inflorescence and details of the spikelet and 
seeds accompany most of the species treatments. The 
photographs are excellent both in quality and 
content, and they supplement the key and the verbal 
descriptions handsomely, clarifying fine points and 
generally improving the likelihood of successful 
identification. They answer a frequent prayer of 
those of us who often find verbal descriptions 
insufficiently descriptive: “oh, please, show me a 
picture of that feature!!” 

The shortcomings I mentioned are exceptions in 
an otherwise splendidly comprehensive treatment of 
the grasses of Oregon and Washington, a book I will 
take with me wherever I go in the Pacific Northwest 
and will supplement my efforts in grass identification 
in California. 

This field guide will be useful long into the future. 
Even after climate change leaves its awful marks on 
the landscape and taxonomists find reason to change 
all the grass names again, the images and descrip- 
tions will connect the changes to what you see in the 
field. While California botanists should consult local 
keys to avoid overlooking species not found in 


O’DONNELL: BOOK REVIEW 113 


Oregon and Washington, this book is worth buying 
for the images and descriptions, serving as an 
invaluable supplement to any California flora or 
field guide. 


—RICHARD O’DONNELL, 1317 Cornell Ave., Berkeley, 
CA 94702. 
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Due to transcription errors, the spatial coordinates 
presented in the Paratypes section are incorrect. A 
revised version of this section is presented below. 

Paratypes (arranged alphabetically by collector 
within county. Note: locality information, including 
spatial coordinates (NAD 83), is provided only when 
given on the collection label.): USA., CALIFOR- 
NIA, Humboldt County: Brush Mountain, 26 May 
1983, Jimerson sn (CAS); FS 5N40 about 3.2 miles 
from NF 60, 40.76532, —123.51479, 4092 ft, 3 Jun 
2013, Mesler 1481] (HSC); FS 6, 1 mile south of Sims 
Mtn, 40.72484, —123.60366, 3600 ft, 5 May 2018, 
Mesler 1665a (HSC); FS 6 near Manzanita Ranch, 
40.71958 —123.54644, 1680 ft, 19 May 2018, Mesler 
1667 (HSC); Friday Ridge Road, near intersection 
with Hwy 299, 7 May 1978, Reed 45 (HSC); Brannan 
Mountain, 4000 ft, 7 July 1911, Tracy 3414 (UC); 
Shasta County: Basin Gulch Camp, 40.35413, 
—122.96176, 2570 ft, 4 May 2013, Mesler 1476 


(HSC); Trinity County: Burnt Ranch USFS camp- 


ground, 412 m, 30 April 1983, Janeway 245 (HSC); 
Weaverville, Kleeberger sn (CAS); Hyampom Rd, 
west of Hayfork, 40.57720, —123.25116, 2100 ft, 7 
June 2011, Mesler 850 (HSC); Wildwood Rd, near 
Gemmill Gulch, 40.43464, —123.05697, 3290 ft, 30 
May 2010, Mesler 818 (HSC);Weaverville, Garden 


Gulch Trail, 40.75234, —122.94912, 2300 ft, 7 June 
2010, Mesler 821 (HSC); North of Hyampom along 
NF 60, 40.57954, —123.25527, 2300 ft, 29 May 2011, 
Mesler 846 (HSC); south of Indian Valley Guard 
Station along NF 7, 40.46557, —123.29167, 4000 ft, 8 
Jun 2011, Mesler $52 (HSC); North of Indian Valley 
Guard Station on NF 10, 40.53945, —123.34832, 4050 
ft, 8 Jun 2011, Mesler 853 (HSC); Philpot Camp near 
Peanut, 40.46500 —123.19028, 2600 ft, 26 May 2013, 
Mesler 1479 (HSC); FS 30NO01 near Wildwood, 
40.40018, —123.00084, 4380 ft, 26 May 2013, Mesler 
1480 (HSC); Ewing Reservoir, Hayfork, 40.56443, 
—123.16867, 2440 ft, 20 May 2018, Mesler 1671 
(HSC); south of Wildwood on FS29N28, 40.36220, 
—123.07039, 4100 ft, 27 Jun 2018, Mesler 1690 
(HSC); FS 60 6.6 mi south of Hwy 299, 3500 ft, 30 
May 1980, Ne/son 5347 (HSC); FS 60 10 mi south of 
Hwy 299, 4400 ft, 30 May 1980, Nelson 5360 (HSC); 
2 miles east of Burnt Ranch, 16 Jun 1964, Spellenberg 
438 (HSC); east side of Underwood Mtn Rd (5N09), 
ca. 4.9 miles south of Hwy 299, 12 June 1980, Taylor 
2855 (CHSC); near summit of Hayfork Mountain on 
Weaverville Rd, 30 June 1923, Tracy 6444 (UC); 
Burnt Ranch, 1500 ft, 27 April 1924, Tracy 6655 
(UC); north of Philpot campground, 2600 ft, 13 May 
1979, York 256 (HSC). 
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